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ABSTRACT 

Wc present Lya luminosity function (LF), clustering measurements, and Lya line profiles based on 
the largest sample, to date, of 207 Lya emitters (LAEs) at z = 6.6 on the 1-deg 2 sky of Subaru/XMM- 
Newton Deep Survey (SXDS) field. Our z — 6.6 Lya LF including cosmic variance estimates yields 
the best-fit Schechter parameters of 0* = 8.5+2° x 10~ 4 Mpc -3 and L\ jya = 4.41°,;® x 10 42 erg s" 1 with 
a fixed a = —1.5, and indicates a decrease from z = 5.7 at the > 90% confidence level. However, 
this decrease is not large, only ~ 30% in Lya luminosity, which is too small to be identified in the 
previous studies. A clustering signal of z = 6.6 LAEs is detected for the first time. We obtain the 
correlation length of ro = 2 — 5 h^ Mpc and bias of b = 3 — 6, and find no significant boost of 
clustering amplitude by reionization at z — 6.6. The average hosting dark halo mass inferred from 
clustering is 10 10 — 10 11 Mq, and duty cycle of LAE population is roughly ~ 1% albeit with large 
uncertainties. The average of our high-quality Keck/DEIMOS spectra shows an FWHM velocity 
width of 251 ± 16km s" 1 . We find no large evolution of Lya line profile from z = 5.7 to 6.6, and 
no anti-correlation between Lya luminosity and line width at z = 6.6. The combination of various 
reionization models and our observational results about the LF, clustering, and line profile indicates 
that there would exist a small decrease of IGM's Lya transmission owing to reionization, but that 
the hydrogen IGM is not highly neutral at z = 6.6. Our neutral-hydrogen fraction constraint implies 
that the major reionization process took place at z > 7. 

Subject headings: galaxies: formation — galaxies: high-redshift — galaxies: luminosity function — 
cosmology: observations 



1. INTRODUCTION 

Understanding physical process of cosmic reionization 
is one of the major goals in astronomy today. Although 
the increase of Gunn-Peterson (GP) optical depths ma y 
be contiguous from low z to z > 6 (jBecker et alJ :2007). 
the evolution of GP optical de pths clearly show the steep 
rise at z ~ 6 towards high-z dFan et al.l 120061 ). On the 
other hand, iDunklev et al.l (|2009f ) find that the polariza- 
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tion data of WMAP place the constraints that instanta- 
neous reionization at the late epoch below z = 8.2 (6.7) 
is rejected at the 2cr (3c) level, and claim that the reion- 
ization process wou ld be extended at z ~ 6 — 11 (see also 
lLarson et af]l2010l for the latest WMAP7 results). How- 
ever, physical models would not easily reproduce such a 
long extended reioniza tion due to the rapid recombina- 
tion of hydrogen (e.g. iFukugita fc Kawa saki 1994; Cen 
120031 ). Observational measurements on neutral hydrogen 
fraction of inter-galactic medium (IGM) at z ~ 6 — 11 
are the missing pieces in this cosmological puzzle. 

Studies of galaxies near the epoch of reionization 
(EoR) at z > 6 are essential not only for understand- 
ing cosmic reionization process but also galaxy forma- 
tion history. The combination of Subaru and VLT 
wide-field cameras and the newly-installed HST/WFC3 
has identified a definitive decrease of UV-continuum 

luminosity function (LF) from z = 6 to 7 — 8 

(lOuchi et all l2009bt lOesch et all 120101 iBouwens et all 



2010al iMcLure et all 120101; ICastellano et al 



Wilkins et alJ|2010b iBun ker et al.ll2009t lYan et all 12009: 



10: 
10: 



2010; 



Wilkins et all l2010t see also Iffickev et " al.l 120091 ). which 
are also reproduced b y recent hydrodynamic simulations 
(|Finlator et al.l l2010f ). Accordingly the star- formation 
rates drop roughly by an order of magnitude from its 
peak of z = 2 — 3 to z = 7, and this decrease implies 
that observations are touching initial formation epoch 
of galaxies. Because galaxies are thought to be sources 
of reionization, the decrease of UV LF would indicate 
that galaxies produce less UV ionizing photons towards 
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TABLE 1 

Imaging Observations and Data 



Band 



Field Namc(s) 



Exposure Time 
(sec) 



PSF size 3 
(arcsec) 



Area 
(arcmin 2 ) 



m-lim 
(3<7 AB mag) 



Date of Observations 



^£921 
^£921 
^£921 
^£921 
^£921 



SXDS-C 
SXDS-N 
SXDS-S 
SXDS-E 
SXDS-W 



30000 
37800 
37138 
29400 
28101 



0.73 (0.81) 
0.75 (0.85) 
0.79 (0.83) 
0.65 (0.83) 
0.75 (0.83) 



590 
734 
783 
610 
521 



26.2 2005 Oct 29, Nov 1, 2007 Oct 11-12 

26.4 2005 Oct 30-31, Nov 1, 2006 Novl8, 2007 Octll-12 

26.2 2005 Aug29, Oct29-30, Novl,2006 Novl8,2007 Octl2 

26.2 2005 Oct 31, Nov 1, 2006 Nov 18, 2007 Oct 11-12 

26.2 2006 Nov 18, 2007 Oct 11-12 



Archival broad-band data c 



B SXDS-C,N,S,E,W 19800 - 20700 0.78 - 0.84 915 - 979 28.1 - 28.4 

V SXDS-C,N,S,E,W 17460 - 19260 0.72 - 0.82 915 - 979 27.7 - 27.8 

R SXDS-C,N,S,E,W 13920 - 14880 0.74 - 0.82 915 - 979 27.5 - 27.7 

i' SXDS-C,N,S,E,W 18540 - 38820 0.68 - 0.82 915 - 979 27.5 - 27.7 

z' SXDS-C,N,S,E,W 11040 - 18660 0.70 - 0.76 915 - 979 26.4 - 26.6 



a FWHM of PSFs in the reduced image. Values in parenthesis indicate the FWHMs of PSF that are matched with the broad-band images in each field. 

Limiting magnitude defined by a 3cr sky noise in a 2"-diametcr circular aperture. 
c The archival broad-band data of SXDS presented in lFurusawa et al.l ( 120081) ■ We show the properties of the 5-subficld images on a single line. Note that 
the exposure time is not a total of the 5 subfields, but 1 subficld, i.e. integration per pixel. More details arc presented in Table 2 of lFurusawa et all J2008T) . 



high redshifts. The production rate of UV ionizing pho- 
tons is close to balance with the recombination rate of 
hydrogen IGM at z ~ 7 — 8, and the reionizing epoch 
may be near these redshifts. The other interpretation of 
the decrease of UV LF is that a moderately high ion- 
izing photon escape fraction, / csc > 0.2, is required to 
keep the universe ionized at z ~ 7 (e.g. lOuchi et al.1 
12009b! : iBunkeret all 120091 ) . The faintest HST/WFC3 
sources at z ~ 7 — 8 show a very blue UV continuum 
slope possibly consistent with extremely young, metal- 
poor stellar populations, whic h is also suggestive of l arge 
esca pe fraction of ,f esc > 0- 3 (jBouwens et alj|2~010bt see 
also iFinkelstein et al.l 120 091) . although there are claims 
that their UV continuum slope measurements include 

otentially large statistical a nd systematic uncertainties 

Schaerer fc de Barrosll2010( ). 

On the other hand, neutral hydrogen of IGM absorbs 
Lya emission line from galaxies via Lya damping wing, 
and dim Lya luminosity, which w ould be identified in 
the evolution of Lya emitter s (e.g. [M alhotra fc Rhoads l 
2004: Kashikawa ct al. 2006; Iv e et al.l 120061: lOta et"afl 
2001) IMalhotra fc Rhoadsl (pOOlT and IHu fc Cowid 
(|2006[ ) find that no significant change of Lya luminos- 
ity function (LF) from z = 5.7 to 6.6, while the study of 
Kashikaw a et al.l (|2006l ) claims that Lya LF evolves from 
z — 5.7 to 6.6. Similarly, Lya LF evolution from z — 5.7 
toward higher redshifts, z ~ 7—7.7, is also under debate. 
Lya LFs at these redshifts are estimated with a reliable 
but o nly 1-2 LAEs at z = 7.0 (jive et al.ll20"ol lOta et all 
l2008f ) and with relatively less-reliable z ~ 7.7 candidates 
with no spectros copic confirmation (jHibon et al.l 120101 : 
iTilvi et ail |2010|) . Moreover, beyond this epoch, only 
weak upper limits are placed on Ly a LF at z ~ 8.8 
(IWillis fc Courbinl [20051 iCubv et al.l I2007J: IWillis et al.l 
I2008t ISobral et all 12001 cf. IStark et al.1 120071b There 
are two conclusions from these z = 6.6 — 7.7 studies, 
no evolution and a decrease of Lya LF, which do not 
agree with each other. These different conclusions of 
Lya LF evolution may be raised by contamination, small 
statistics and systematic errors such as cosmic variance. 
Moreover, a careful argument is needed for the interpre- 
tation of Lya LF evolution, because properties of star- 
formation galaxies are changing, which are already found 
in the decrease of UV LF of dropout galaxies as dis- 



cussed above. It is also important to constrain reion- 
ization and galaxy formation models with other observa- 
tional quantities. The stacked spectra of z = 6.6 LAEs 
show no clear signal of Lya damping wi ng absorption 
(jKashikawa et al.ll2006l ; IHu fc Cowid l2006f ). Recent the- 
oretical studies predict that clustering measurements of 
LAEs can be an independent probe of reionization. Since 
Lya lines of galaxies residing in ionized bubbles selec- 
tively escape from a partially neutral universe, cluster- 
ing ampli tude of observed LAEs would be boos t ed at 
the EoR jlFurlanetto et all 120061: iMcQuinn et al.l 120071 : 
iLidz et al.ll2009l ; cf. Illiev et al.ll2008f ). Although the im- 
portance of z > 6 LAE clustering is claimed from the 
theoretical studies, no ob servational study has pr ovided 
a reliable measurement. IKashikaw a et al.l ([2006D mea- 
sured angular correlation function (ACF) of 58 LAEs at 
z = 6.6 which is so far the largest sample of z = 6.6 
LAEs, but no significa nt signal is detected, p robably due 
to the small statistics (jMcQuinn et al.ll2007t) . 

To address these issues in statistics of z = 6.6 LAEs, 
we are conducting an extensive survey that is the largest 
ever perfor med in terms of area and number of ob- 
jects (cf. IMalhotra fc Rhoadsl 12004 IHu et all 120051 ; 
IKashikawa et al.l |2006[ K Our survey field is a large 1 
deg 2 area of Subaru/XMM- Newton Deep Survey (SXDS) 
fie ld well-separated from t he Subaru Deep Field (SDF) 
of IKashikawa et all (|2006l ). The combination of SXDS 
and SDF data will allow us to provide a large statis- 
tical sample with cosmic variance error estimates. The 
spectroscopic component of this survey is underway with 
Keck/DEIMOS and Magellan/IMACS (Ouchi et al. in 
prep.). The aim of this study is to supply the reliable 
statistical measurements of z = 6.6 LAEs based on our 
large sample that are useful for constraining the existing 
and forthcoming cosmological models of cosmic reion- 
ization and galaxy formation. Our z = 6.6 LAE re- 
sults from the large statistics will also be good baselines 
for on-going and future studies of LAEs at z > 7 per- 
formed with the ex isting nea r-infrared imagers o f VLT 
(IWillis et al.1 120081). S ubaru (jTokoku et al.ll2008D . and 



VISTA (|Nilsson et al.|[2007D and the future facilities of 
JWST, E-ELT, GMT, and TMT. We present our survey 
and sample in 52 luminosity functions in Ej3j cluster- 
ing properties in 21 Lya line profiles in Sj5l and discuss 
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Fig. 1. — Color magnitude diagram of narrow-band excess color 
[z! - iV£921) vs. narrow-band magnitude (JVB921). Black dots 
present colors of all the detected objects. Black filled circles and 
triangles denote our z = 6.6 LAEs. Specifically, the triangles in- 
dicate LAEs with z'-band magnitudes fainter than that of the 3cr 
level, and show their \a lower limits of their z' — NB921 colors. 
Red and blue open symbols mark spectroscopically-identified ob- 
jects in the redshift range of LAEs and interlopers, respectively. 
We define the redshift ranges of LAEs as 6.45 — 6.65 (see Figure 
[4] for the redshift range) . The green lines indicate 2a errors of the 
color of z' - NB921 for a source with a color of z' - NB921 = 0.04, 
which corresponds to the median color of all objects. Dashed and 
dotted lines represent our color cut for narrow-band excess and the 
1<t limit of z' data, respectively. Note that 7VB921 magnitudes are 
total magnitudes, while the colors of z' — NB921 are defined with 
a 2"-diameter aperture. 

reionization and galaxy formation by comparisons with 
theoretical models in SJBJ Throughout this paper, magni- 
tudes are in the AB system. We adopt a cosmology pa- 
rameter set of (h, Q m , flA,n s ,as) = (0.7, 0.3,0.7, 1.0,0.8 
consistent with the WMAP results (|Komatsu et alj|200~ 
lLarson et aD2010t ). 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Imaging Observations 

We carried out extensive deep narrow-band ima ging 
with Subaru/Suprime-Cam (iMivazaki et al.l I2002T ) for 
SXDS in 2005-2007. We summarize details of our ob- 
servations as well as image qualities in Table [T] We 
used the narrow-band filter, NB921, with a central wave- 
length of A c = 9196A and a FWHM of 132A. The re- 
sponse curves of the A^S921-band and broad-band fi l- 
tcrs are presented in Figure 1 of lOuchi et al.l (|2008[) . 
Five pointings of Suprime-Cam covered 5 subfields o f 
SXDS, SXDS-C, N, S, E, and W dFurusawa et al.ll2008l) . 
by 8.3, 10.5, 10.3, 8.2, and 7.8 hour on-source integra- 
tion, respectively. In addition to these narrow-band im- 
ages, we use archival data of deep broad -band (B, V, R 
i! an d z') images of the SXDS project (jFurusawa et al.1 
l2008h as summarized at the second half of Table [1] The 
narrow-band data are reduced with the Su prime-Cam 
Deep field REDuction p ackage (SDFRED; lYagi et al.l 
I2002t lOuchi et alll2004ah . With the standard parame- 
ter sets of SDFRED, we perform bias subtraction, flat- 



FlG. 2. — Two color diagram of narrow-band excess color 
(z' — iVi?921) vs. continuum color («' — z'). Black dots present 
colors of all the detected objects. Our z = 6.6 LAEs are shown by 
black and red filled circles+triangles. The triangles have a z'-band 
magnitude fainter than 3<r level, and present the la lower limit 
of z' — NB921. Red and blue circles+triangles indicate spectro- 
scopically identified objects in and out of the LAE redshift range, 
respectively. We take the same LAE redshift range as that of Fig- 
ure [T] For display purpose, objects with no color measurements 
(no z'-band detection at the 3a level) are placed at a color of 
i' — z' = 1.3. The vertical and horizontal solid lines show our 
color criteria (eq. [T) for sources with a detection in z' band. 
Colors are defined with a 2"-diameter aperture. Curves present 
tracks of model interloper galaxies at differe nt redshifts. Green 
lines show 6 templates of starburst galaxies ((Kinney et af] |1996l ) 
up to z = 2, which are 6 classes of starburst galaxies with dif- 
ferent dust extinction (E(B — V) = 0.0 — 0.7). The narrow-band 
excess peaks in the green lines correspond to the emission lines 
of Ha(z = 0.4), [Om](z = 0.8), H/3(z = 0.9), or [On] (2 = 1.5). 
Black lines represent colors of typical elliptica l, spiral, and irregu- 
lar galaxies (Coleman. Wu, & Wccdman 1980) which arc redshifted 
fro m z = to z = 2. Yellow star marks are 175 Galactic stars given 
bv lGunn fc Strvkerl (IT983T) . 

fielding, distortion+atmospheric-dispcrsion corrections, 
sky subtraction, image alignments, and stacking. Before 
stacking, we mask out bad data areas such as dead pixels 
and satellite trails. Cosmic rays are removed in the pro- 
cess of stacking with the rejected-mean algorithm. The 
final NB921 images have the seeing size of 0".7 — 0".8 
in FWHM, and reach the 3<7 limiting magnitudes of 
26.2 — 26.4 mag in a 2".0-diameter aperture. 

We use neither contaminated areas with halos of bright 
stars and CCD blooming nor low signal-to-noise (S/N) 
ratio regions located around the edge of the FoV, which 
are caused by dithering. After we remove these low- 
quality regions from our catalog f £|2.2[) . the effective to- 
tal areas are 3238 arcmin 2 . This effective area cor- 
responds to the survey volume of 8.0 x 10 5 Mpc 3 at 
z = 6.565 ± 0.054, if we assume a simple top-hat se- 
lection function of LAEs whose redshift distribution is 
defined by the FWHM of our narrow-band filter. 

During the observations, we took images of spectropho- 
tometric standard star of GD71 with iVB921-band filter 
(|Bohlin et al.lll995l ) . The standard star was observed un- 
der photometric condition in 2006 November 18 and 2007 



Ouchi ct al. 



October 11-12. We calculate photometric zero-points 
from the standard star data. We check these photo- 
metric zero points based on colors of stellar objects in 
our field and 175 Galacti c stars calculated from spectra 
given in lGunn &: Strvken (J1983I ). We find that colors of 
stellar objects in our da ta are consistent with those of 
iGunn fc Strvkert s (|1983| ) stars within ~ 0.05 magnitude. 
Our narrow-band images are aligned with the deep op- 
tical BVRiz images from the SXDS survey (Table QJ 
iFurusawa et al.l [20081) based on hundreds of stellar ob- 
jects commonly detected in both narrow- and broad-band 
images. The astrometry of our objects is the same a s 
those of SXDS version 1.0 catalog (jFurusawa et al.ll2008l ). 
The errors in the relative positions of objects are ~ 0".04 
in r.m.s. Th e r.m.s accuracy of the absolute positions is 
estimated in IFurusawa et all (|2008D to be ~ 0".2. After 
the image registration, we homogenize the PSF sizes of 
broad and narrow-band images, referring to these stel- 
lar objects. The PSF sizes of our narrow-band images 
match to those of broad-band images with an accuracy 
of AFWHM ~ 0".01. 

2.2. Photometric Sample of z = 6.6 LAEs 

We identify 286,510 objects in the N B921 images down 
to_/V\B921 = 26.0 with SExtractor (jBertin fc Arnoutsl 
[19961) . We measure both MAG_AUTO of SExtrac- 
tor and 2".0-diameter aperture magnitudes. We adopt 
MAG_AUTO as total magnitudes, while we use a 2".0- 
diameter aperture magnitude to measure colors of ob- 
jects in order to obtain colors of faint objects with a 
good signal-to-noise ratio. We correct the magnitudes 
of objects for Galactic extinction of E(B — V) = 0.020 
([Schleeel. Finkbeiner. fc Davislll998l ). 

We plot a color-magnitude diagram in Figure [1] 
for our objects. Figure Q] shows narrow-band ex- 
cess color, z' — NB921, and narrow-band magnitudes, 
NB921. Figure [2] presents a two-color diagram based 
on the iV_B921-detcction catalog, together with colors 
of model galaxies and Galactic stars. We plot col- 
ors of 3,249 spectroscopically-identified objects which 
include our 16 z = 6.6 LAEs (s ee 3Q) and 3,23 3 
foreground/background interlopers ([Ouchi et al.l 120081 ). 
Spectroscopically-identified LAEs are located in the 
upper-right part of the two-color diagram, having a 
narrow-band excess of Lya emission (z! — NB921) and 
red continuum colors of GP trough (i' — z'). 

Based on the color diagram, we select candidate LAEs 
with the narrow-band excess, no detection of blue con- 
tinuum flux, and the existence of GP trough, by the color 
criteria, 

z' - NB921 > 1.0 and B > B 2 <r and V > Vi a and 

[(*' < z' Za and i>-z>> 1.3) or (z' > z' 3<7 )} , (1) 

which are similar t o those in the study of SDF 
([Taniguchi et al.ll2005l: iKashikawa et alJl2006( ). B 2a and 
Via arc defined as 2a limiting magnitudes of B and V 
bands, respectively (B^a = 28.7 and Vi a = 28.2), which 
ensures no detection of continuum bluer than Lyman 
break (~ 6900A) for objects at z = 6.6 
detection limit (z' 3a = 26.5). 

We obtain a photometric sample of 207 LAEs at z ~ 
6.56 ± 0.05 down to JVJ392 1 = 26.0 in a comoy ing survey 
volume of 8 x 10 5 Mpc 3 (|Ouchi et al.ll2009al) . Our se- 
lection criteria correspond to LAEs with the rest-frame 
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Fig. 3. — Spectra and snapshots of our z = 6.6 LAEs confirmed 
with Keck/DEIMOS. Each object has a spectrum in the right panel 
and snapshots of B, V, R, i', z' , and NB921 images in the left 
panels. Each snapshot is presented in a 6" X 6" box. The object 
name and redshift are presented in the left and right corners of 
each spectrum panel, respectively. The right bottom panel shows 
a typical DEIMOS spectrum of the sky background that is obtained 
in the process of sky subtraction. 

equivalent width, EWq, of > 36A, if a flat continuum 
spectrum (/„ =const) is assumed. If the realistic spec- 
trum of LAEs with a Gunn-Peterson trough is assumed, 
the limit of rest-frame equivalent width is EWq > 14 A 
and a limiting line flux of / > 5 x 10~ 18 erg s _1 cm~ 2 
corresponding to L > 2.5 x 10 42 erg s _1 at z ~ 6.56. 

2.3. Spectroscopic Observations 

We conducted deep spectroscopic follow-up for our 
z = 6.6 LAEs with Keck/DEIMOS (|Faber et al.|[200l 
on 2007 November 6, 2008 August 1, and 2008 Octo- 
ber 3. We took spectra of 3.0 hour on-source integra- 
tion for 3 masks and 2. 3- hour integration for 1 mask 
with a 1".0 slit and the 830G grating, which covered 
~ 6000 — 10000 A and gave a medium high spectral res- 
olution of R ~ 3600 at 9200A. We observed 30 out of 
207 LAEs including faint LAEs whose expected flux is 
below our observational limit, and obtained 19 spectra 
with significant signals. All of the 19 spectra show a sin- 
gle line at around 9200A with no detectable continuum. 
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TABLE 2 
Lyo Emitters with spectroscopic redshifts 



Object Name 
(1) 


a(J2000) 

(2) 


<5(J2000) 
(3) 


Z 

(4) 


m NB 

(mag) 

(5) 


(mag) 
(6) 


L(Lya) 

(10 43 erg s" 1 ) 

(7) 


AVpwHM 

(km s" 1 ) 
(8) 


EW 

(A) 

(9) 


M uv 

(mag) 

(10) 


SFR 

(Moyr- 1 ) 
(11) 


NB921-C-106098 
NB921-C-22057 
NB921-C-34609 


2:17:03.4908 
2:18:20.6670 
2:18:19.3797 


-4:56:19.158 
-5:11:09.664 
-5:09:00.550 


6.589 
6.575 
6.563 


24.7 
24.9 
24.8 


26.8 

> 27.7 

27.3 


1.2 ±0.1 

1.0 ±0.1 
0.9 ±0.2 


269 ± 37 
304 ± 58 
247 ± 98 


105.9^ 5 .4 
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~ -20.2 

< -19.9 
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-21.1 ±0.7 


15.0 ±9.7 


NB921-C-43803 


2:18:26.9734 


-5:07:26.828 


6.554 


25.0 


26.9 


0.7 ±0.2 


199 ± 190 


on ry-\-o2o . 2. 

oz.z_ 21 3 


~ -20.7 


~ 10.3 


NB921-C-50823 a 
NB921-C-65683 
NB921-N-71598 
NB921-N-77765 


2:17:02.5772 
2:17:01.4777 
2:18:44.6534 
2:18:23.5317 


-5:06:04.480 
-5:03:09.291 
-4:36:36.501 
-4:35:24.193 


6.545 
6.493 
6.621 
6.519 


24.8 
24.7 
24.7 

24.7 


26.9 
26.1 
27.4 
26.4 


0.9 ±0.3 
3.5 ±12.7 
2.6 ±0.3 
0.4 ±0.8 


257 ± 69 
329 ± 70 
367 ± 46 
266 ± 80 


52.9_ 35 g 

^ ioo 
Z ioo 

45 2+°° 
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~ -20.5 

< -19.9 
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<5 

<5 

22.9 ±14.8 


NB921-N-79144 
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24.1 
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0.9 ±1.2 
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42.7 ±15.7 


NB921-W-25755 
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-4:55:56.703 
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25.1 


> 27.7 


0.7 ±0.2 


134 ± 42 


l«-llu6.8 


< -19.9 


<5 


NB921-W-30717 
NB921-W-30999 a 
NB921-W-31790 a 
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2:16:53.9169 


-4:55:56.788 
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-5:06:01.379 


6.617 
6.512 
6.590 


25.3 

25.1 
25.8 
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1.3 ±0.2 
1.3 ±0.8 
0.4 ±0.1 


246 ± 44 
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142 ± 154 


153.01-1.4 
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45 9+°° 


~ -20.1 

< -19.9 

< -19.9 


~5.9 
<5 
^5 


NB921-W-35578 


2:16:50.9037 


-5:10:16.221 


6.535 


25.4 


> 27.7 


0.6 ±0.3 


279 ± 210 


104.4l-. 


< -19.9 


$5 


NB921-W-36217 


2:16:50.4202 


-5:05:45.130 


6.586 


25.5 


27.3 


0.5 ±0.2 


128 ± 145 


34.7±^ 3 5 ' 6 


< -19.9 
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Because we allowed redundant observations for 6 LAEs 
that can be included in two different masks. Three out 
of 19 detected objects were taken for the same targets. 
Thus, totals of the observed and identified LAEs are 24 
(= 30 — 6) and 16 (= 19 — 3), respectively. We hereafter 
refer to the 24 galaxies as the spectroscopic sample. The 
3 duplicate spectra with the signals present an emission 
line at the same wavelength as the original spectra, and 
confirm that the emission lines are real signals. We use 
these 3 duplicate spectra for the stacking analysis in Sj5] 
We present our spectra in Figure [3J We have con- 
firmed that no spectra show signatures of an [Oiii] 5007 
emission line (at ~ 7000A) and an [On] 3727 emis- 
sion line (at ~ 5200A) from a z = 0.40 Ha emitter or 
an [On] emission line (at ~ 6800A) from a z = 0.84 
[Oiii] emitter, and found that these objects are nei- 
ther a foreground Ha nor [Oiii] emitters. If there is 
a z ~ 1.47 [On] emitter, it should show a AA3726,3729 
doublet whose separation is 7A. Our DEIMOS spectra 
with a FWHM spectral resolution of 2.6A can distin- 
guish the doublet from a Lya line with a characteristic 
asymmetric profile. Our DEIMOS spectra confirm no 
such signature of [Oil] doublet, but mostly clear asym- 
metric line profiles with an extended red wing that is 
typical for a high-z Lya line. One exception is the spec- 
trum of NB921-W-35578 that has three spikes in the 
range of ~ 9156 — 9162A. The two-dimensional spectrum 
of NB921-W-35578 presents a significant residual of sky 
subtraction around ~ 9159 — 9164A, although a signal 
of emission line is clearly seen in the wavelength range 
of free of residuals. The complicated spectral shape of 
NB921-W-35578 is probably made by errors of sky sub- 



traction. Table [5] summarizes properties of the 16 LAEs 
in our spectroscopic sample. The Lya emission lines 
of these LAEs are strong, EWq > 20, which are con- 
sistent with the expected EWq from our color criteria 
[EW > 14A; 32~2"j) . Figure H plots the redshift distribu- 
tion of our spectroscopic sample, together with a selec- 
tion function simply estimated from the filter response 
curve of 7V.B921 band. The estimated selection function 
is similar to the redshift distribution of our spectroscop- 
ically identified LAEs. 

3. LUMINOSITY FUNCTION 

3.1. Completeness and Contamination of our Sample 

We estimate completeness and contamination of our 
sample in the same manner as lOuchi et al.l (|2008l ). First, 
we estimate the detection completeness of narrow-band 
images, /det, as a function of narrow-band magnitude. 
We distribute 7400 artificial objects with a PSF profile 
that mimic LAEs on our original 3238 arcmin 2 images 
after adding photon noise, and detect them in the same 
manner as for the detection of our LAE catalogs with 
SExtractor. We repeat this process 10 times, and com- 
pute the detection completeness. We plot the detection 
completeness as a function of magnitude in the top panel 
of Figure [5j The top panel of Figure [S] presents that the 
detection completeness is typically ~ 80 — 90% for rel- 
atively luminous sources with A.B921 < 25.5. The de- 
tection completeness is ~ 50% in the faintest magnitude 
bin of our sample, NB921 = 25.5 - 26.0. 

Second, we estimate the contamination of our LAE 
samples. There are no interlopers in our spectroscopic 
sample. We also use a spectroscopic cata log of 3,233 
SXDS objects at low redshifts of z = — 6 (|Ouchi et al.l 



Ouchi ct al. 



10 



^ 







6.45 



6.55 



6.65 



spec 



Fig. 4. — Redshift distribution of our LAEs with a spectro- 
scopic identification. Histogram presents LAEs confirmed by our 
Keck/DEIMOS observations. Dashed line represents the selection 
function of LAEs that is simply calculated from the response curves 
of NB921 filter. The selection function is normalized by the num- 
ber of the identified LAEs. 



l2008f ). and find none of these 3,233 low-z objects are 
included in our z = 6.6 LAE sample. We define the 
contamination fraction, / con t> with 



/cont = -/Vlowz/iVall 



(2) 



where Mowz and N a u are the numbers of low-z interlopers 
and all objects, respectively, in our spectroscopy sample. 
The totals of the spectroscopically observed and iden- 
tified LAEs are 24 and 16, respectively f §2.3|) . Since 
there are 8 objects with no identification in our spectro- 
scopic sample (see Section |2~3|) . we calculate / cont for the 
following two extreme cases. If all of these unidentified 
objects are real LAEs whose Lya lines are simply too 
faint or extended to be detected in our spectroscopy, we 
find Ni owz /N a ii = 0/24 because of no interlopers in our 
spectroscopic sample. If all the unidentified objects are 
interlopers, N\ owz /N a ii = [24— 16]/24. Thus, the contam- 
ination rate is taken within the range of f con t ~ — 30% 
for our LAE samples. Note that there are no obvious con- 
taminants in our follow-up spectroscopy and the SXDS 
catalogs. Moreover, the on-going DEIMOS and IMACS 
spectroscopy have already confirmed 70% of our LAEs 
at log-Liyya > 43.0 erg s , and found no contaminants 
in our sample (Ouchi et al. in preparation). Although 
the effect of contaminants may be negligible, we use this 
upper limit of contamination (/ con t — 30%) to derive the 
upper limits of our measurements in <gj 

3.2. Surface Number Density and Cosmic Variance 

The bottom panel of Figure[5]shows surface densities of 
LAEs and all narrow-band detected objects (designated 
as All Objects'). Red circles are the average surface 
densities. Black points with 5 different symbols indicate 
the surface densities in the 5 subfields (~ 0.2 deg 2 ) of 
Suprime-Cam, i.e., SXDS-C, -N, -S, -E, -W. The detec- 
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Fig. 5. — Top panel: Detection completeness of our NB921 
images in percentage. Circles, hexagons, triangles, squares, and 
pentagons represent the completeness of a magnitude bin (Am = 
0.5 mag) in the 5 subfields, SXDS-C, SXDS-N, SXDS-S, SXDS-E, 
and SXDS-W, respectively. For presentation purpose, we slightly 
shift all the points along the abscissa. Bottom panel: Surface 
densities of objects detected in the NB921 data. The lower and 
upper points indicate surface densities of our z = 6.6 LAEs and 
all the objects detected in the narrow band, respectively. Black 
circles, hexagons, triangles, squares, and pentagons plot the surface 
densities in the 5 subfields of SXDS-C, SXDS-N, SXDS-S, SXDS- 
E, and SXDS-W, respectively. We distinguish between the raw 
and completeness-corrected surface densities with the open and 
filled symbols, respectively. The red filled circles represent the 
surface density averaged over our entire survey field. The errors 
are given by Poisson statistics for black symbols, while the errors of 
red symbols are the geometric mean of Poisson errors and cosmic 
variances calculated from eq. ((3j . To avoid overlaps of points, we 
slightly shift all the black points along the abscissa with respect 
to the corresponding red filled circles. Right vertical axis indicates 
the numbers of objects, i.e. 7V/(0.5mag), identified in our entire 
survey area. Blue stars denote the surface density of z ~ 6.6 LAEs 
obtained bv lTaniguchi et al.l d2005l ). 

tion completeness correction is applied based on the sim- 
ulation results described in Section I3.ll After the com- 
pleteness correction, surface densities of z = 6.6 LAEs 
among the 5 subfields differ up to by a factor of ~ 2 in 
faint magnitude bins (NB921 = 25 — 26) and a factor of 
~ 10 in bright magnitude bins (NB921 < 25), while the 
difference for All Objects' is negligibly small. Accord- 
ingly, the large variance of z = 6.6 LAEs is not artifacts, 
but real. These large differences of LAE surface densities 
probably come from the combination of cosmic variance 
and Poisson error s. Following the procedures used in 
lOuchi et al.l (|2008h . we evaluate the cosmic variance in 



our survey area, a„ 



with 



2 
Ogtlsf 



Cg — Ogilsf (cDM/o"DM:lsf ) 

((£ g:lsf -S g ) 2 )-S, 



\m 



(3) 

(4) 



Statistics of Lya Emitters at z ~ 7 



o 
ft 



II 

fc£ 
O 

5 



10" 



10" 



10" 



10" 




z = 6.6 



This work (all) 
)ADOThis work (subfield) 
•k Kashikawa ct al. (2006) \ 

Hu et al. (2005) 
Malhotra & Rhoads , (2004) 



42.0 43.0 

logl(Lya) [ergs s" 1 ] 



44.0 



Fig. 6. — Luminosity functions (LFs) of LAEs at z = 6.6, to- 
gether with those of previous measurements for z = 6.5 — 6.6 LAEs. 
Filled red circles denote the best-estimate LF at z = 6.6 from our 
SXDS data. The error bars include uncertainties from statistics 
and cosmic variance (eq. \3$. Black open circles, hexagons, trian- 
gles, squares, and pentagons represent the LFs of five ~ 0.2 deg 2 
subficlds, SXDS-C, SXDS-N, SXDS-S, SXDS-E, and SXDS-W, re- 
spectively. In order to avoid the overlaps of points, we slightly 
shift all the open symbols along the abscissa with respect to the 
corresponding red filled circles. Bl ue star marks present z ~ 6.6 LF 
derived in the ~ 0.2 deg 2 SDF bv lKashikawa et al.l 120061) (photo- 
metric sample), while blue cro sses are th e meas urement and upper 
limit of 2 = 6.5 LF given by IHu et al.l (120051 s ). Blue dotted line 
shows the LF estimated bv lMalhotra fc Rhoads (2004}. Red solid 
line is the best-fit Schechter function to the combination of SXDS 
and SDF LFs at z = 6.6. 

where ctdm and CDM:isf are the rms fluctuation of 
dark matter in all the survey volume and the vol- 
ume of 5 subficlds (~ 0.2 deg 2 ), respectively. We 
calculate the fluctuations of dark matter with the 
po wer spectrum, adopting the transfer fu nction given 
bv lBardeen. Bond. Kaiser, fc Szalavl (|1986lh a g: i SI is the 
fluctuation of number density of LAEs for one subfield. 
S g: isf and Eg are LAE's surface densities in a subfield 
and the entire survey areas, respectively. 

Since these estimates of cosmic variance are based on a 
large but single contiguous field, it is important to check 
whether our field is located at the sky of an overdense or 
underdense region. Moreover, a large-scale overdensity 
or undcrdensity of Lya sources could also be produced 
by an inhomogeneous distribution of Lya absorbers (i.e. 
neutral hydrogen) along the line of sight. In Figure [5] 
we compare the surface densities of our LAEs with those 
selected with the same NB921 fil ter in a completely in- 
depend ent sky of ~ 0.2 deg 2 SD F (jTaniguchi et al.ll2005t 
see also lKashikawa et al.ll2006[ ). We find that the surface 
de nsities of our z = 6.6 L AEs are consistent with those 
of iTaniguchi et all (J2005I ) within the scatters and Pois- 
son errors of our 5 su bfields. Moreover, Figure [5] shows 
that ITaniguchi et alj 's measurements scatter around our 
average surface-density curve. Thus, we conclude that 
our field has neither signature of overdensity nor undcr- 
density, and that our 5 subfields appear to represent the 
average cosmic variance. 



3.3. Lya Luminosity Function 

We deri ve Lya LF o f LAE s at z = 6.6 in the same 
manner as lOuchi et al.l ()2008f ) . We obtain the number 
densities of LAEs in each magnitude bin by dividing the 
observed number counts of LAEs by the effective survey 
volume defined as the FWHM of the NB921 bandpass 
times the area of the survey. Here, we calculate the Lya 
luminosity of each object with the response curves of nar- 
row and broad bands by subtracting the continuum emis- 
sion measured from the continuum magnitude from the 
total luminosity in the narrow band. In this calculation, 
we use total magnitude of narrow-band images. The con- 
tinuum emission is estimated by the narrow-band excess 
color defined with a 2" aperture, so as to keep high signal- 
to-noise ratios and to avoid the effects of source con- 
fusion on broad-band images with high source density. 
Lya luminosity estimates of this photomet ry technique 
i s test ed for z = 3 — 6 LAEs in Figure 15 of lOuchi et al.l 
( 2008) which verifies that the ratios of these luminosity 
estimates to independent spectroscopy measurements are 
about unity on average. 

Although this procedure of Lya LF derivation is a 
simple classical met hod that was taken by most of pre- 
vious studies (e.g. lOuchi et all 120031: lAiiki et~al] l2003t 
IHu et al.lf200l iMalhotra fc RhoadsH2004f ). there are two 
possible sources of uncertainties in Lya LFs derived with 
this classical method. (1) The narrow-band magnitude of 
LAEs of a fixed Lya luminosity varies largely as a func- 
tion of redshift. Thus Lya luminosity may be over- or 
under-estimated for some LAEs. (2) The selection func- 
tion of LAEs in terms of EW also changes with redshift; 
the minimum EW value corresponding to a given (fixed) 
narrow-band excess, z' — NB921, becomes larger when 
the redshift of the object is far from the redshift cor- 
responding to the center of NB921 bandpa ss. To derive 
Lya L F with no bias originated from (l)-(2). IOuchi et all 
(|2008l) have carried out Monte-Carlo simulations with a 
mock catalog of their z = 3 — 6 LAEs uniformly dis- 
tributed in comoving space that are produced by a set 
of the Schechter parameters (a, (f)*,L*) and a Gaussian 
sigm a for a probability distribution of EWq. Figures 16- 
18 of lOuchi et al.l ()2008| ) compare LFs from the unbiased 
simulations with those from the classical method, and 
indicate that the LFs from the classical method are con- 
sistent with those from the simulations. This is because 
the uncertainties (1) and (2) in the classical me thod are 
negligible and/or cancel out (|Ouchi et al.ll2008T) . Thus, 
our simple classical method gives good estimates of Lya 
LF. 

Figure [6] presents the Lya LF of our LAEs at z = 6.6. 
Note that we have corrected the detection-completeness 
by weighting with /dct measured in Section f3.ll To check 
the cosmic variance and the accuracy of our results, we 
plot the estimates of LFs from the entire 3238 arcmin 2 
(red filled circles), together with those from the five sub- 
fields (black open symbols). In Figure El we calculate 
cosmic variances with eq. ((3]), and include these un- 
certainties in the error bars of LFs of the entire field. 
We find that the scatters of the measurements are large 
among the 5 subfields with differences by up to a fac- 
tor of ~ 10, although the typical scatters of the subfield 
results are not far beyond the errors of Poisson statistics. 

We plot previous measurements of Lya LFs at z ~ 
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Fig. 7.— Evolution of Lya LFs up to z = 6.6. Red filled 
circles are the best estimates of z = 6.6 LAEs from the entire 
SXDS sample, and red solid line is the best-fit Schcchtcr function of 
z = 6.6 LAEs. Blue filled circles and solid line are data points and 
th e best-fit Schechter function, respectively, of z = 5.7 LAEs given 
by IQuchi et al.l 1 120081) . Note that the error bars of z = 6.6 and 
5.7 data points (red and blue filled circles) represent uncertainties 
from statistics and cosmic variance. Cy an solid line is the best- 
fit Schechter function of z = 3.1 LAEs dOuchi et al.| [2008). The 
LF decreases from z = 5.7 to 6.6 significantly, while no significant 
evolution can be found between z = 3.1 and 5.7. For comparison, 
we plot LF estimates from each of the five ~ 0.2 deg 2 subfields with 
the same open symbols as found in Figure [6] These open symbols 
illustrate that with the data of a single ~ 0.2 deg 2 field alone (e.g. 
open circles down to logZ/Lya — 42.6) which is a typical survey 
size of previous studies it is difficult to distinguish whether or not 
z = 6.6 LFs show evolution (decrease) with respect to z = 5.7. 
Dashed and dotted lines represent the best-fit Schechter functions 
to our z = 6.6 LF with a <j>* and L* fixed to that of z = 5.7, 
respectively. 

6.6 in Figure [5] The previous measurements of Lya 

L Fs show large scatte rs. If one compares the LFs 

of IMalhotra fc Rhoadsl (120041) fdo _tted lineV lHu et al.l 
(|2005l) (crosses), and iKashikawa et al.l (|2006l) (stars) 
which are derived from a moderate number of LAEs 
and a moderately-wide field, their measurements do not 
appear to agree. This large difference causes a long- 
standing argument of Lya LF evolution between z = 5.7 
and 6.6. However, all of these previous measurements fall 
about within the scatters and errors of LFs derived from 
our 5 subfields each with a ~ 0.2 deg 2 area. Thus, the 
discrepancies of LF measurements between these previ- 
ous studies are probably due to the combination of cos- 
mic variances and small statistics. Because our z = 6.6 
LFs from the 5 subfields cov e r the r ange of measure - 
men ts of IMalhotra fc Rhoadsl (|2004D . iHu et all (|2005l ). 
and IKashikawa et al.l (|2006l ). we conclude that our LF 
measurements agree well with those derived in previous 

studies. 

Schechter function (|Schechterl 119761 ) is fit to the Lya 
LFs composed of our large area data with cosmic variance 
errors and the previous LF estimate in the indepen dent 
876 arcmin 2 area of SDF (jKashikawa et al.ll2006t ). 12 



Fig. 8. — Error ellipses of our Schechter parameters, L? a and 
4>* . Red contours represent the fit of z = 6.6 LF with the fixed 
slope of a = —1.5 based on SXDS and SDF data. The inner and 
outer contours indicate 68% and 90% confidence levels, respec- 
tively, which include cosmic variance erro rs. Blue contours denote 
z = 5.7 LFs given byfOuchi ct al. ( 2 0081 ) , which arc similarly de- 
rived with cosmic variance errors. For a fair comparison with our 
z = 6.6 LF, wc show error ellipses of the 2 = 5.7 LF de rived by the 
classical method (see more details inlOuchi ct al. 2008]). The error 
ellipses of the z = 5.7 LF are larger than those of our z = 6.6 LF. 
This is because the data of z = 5.7 LF have more uncertainties of 
cosmic variance. Black solid and dotted lines indicate 1, 2, and 3 
sigma confidence levels of z = 6.6 and z = 5.7 LFs with no cosmic 
varia nce errors previously d erived solely with the smaller data of 
SDF JKashikawa et al.ll2006T) . 

The Schcchtcr function is defined by 

<t>{L)dL = <j>*(L/L*) a exp{-L/L*)d{L/L*). (5) 

With a total of 265 (= 207+58) LAEs at z ~ 6.6 in a total 
area of 4114 (= 3238 + 876) arcmin 2 data in the SXDS 
and SDF, we obtain the best-fit Schcchtcr parameters of 
(j>* = 8.5±|° x 10- 4 Mpc- 3 and L* hya = 4A±™ X 10 42 
erg s _1 with a fixed a = —1.5, which are summarized 
in Table [3] Because the difference in x 2 for a values is 
insignificant, we fix a to —1. 5, which is a fiducial value 
used for low-z Lya LFs (e.g. IMalhotra fc Rhoadsl 120041 : 
IKashikawa et all 120061 : lOuchi et al. 2008ft . The best-fit 
Schcchtcr function is shown in Figure|H]with the red solid 
line. We also estimate the best-fit Schechter parameters 
with a fixed a = —1.7, because the value of q~ —1.7 is 
recently reported for z = 2— 6 LAEs (|Cassata et al.ll2010l : 
see also iRauch et al.1 120081) as well as z ~ 7 dropouts 
fe.g JOuchi et al.ll2009bHOesch et al.ll2010l : lMcLure et al.l 
120101 ). The best-fit parameters for a = —1.7 are <p* - 



6.9l 2 'j x 10- 



^pc" 3 and L 



Lya 



4.9 



+0.9 



x 10 erg s 



Twenty six out of 29 LAEs used in Malhotra & Rhoads ( 2004) 



Note that our data point at the bright end (log L(Lya) = 
43.5 erg s _1 ) appears to exceed the best-fit Schechter 
function in Figure|Hl although the data point is consistent 

are provided from a su bsamplc of the SDF data IITaniguchi et all 
120051 ; IKashikawa et aI.II2006 | ). To avoid using the same LAE data, 
wc do not include Malhotra & Rhoads ( 2004|) data points for our 
Schechter fitting. We also do not use data points in papers that 
are not published yet in a refereed journal. 
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TABLE 3 

L/Ya LUMINOSITY FUNCTION. 



^Lya 



(10- 4 Mpc- 3 ) (lO^ergs" 1 ) 
(1) (2) (3) 



(1) 



(5) 



(10 



n obs 

- 4 Mpc" 3 ) (10 3 
(6) 



obs 
l J L,ya 

'erg s — 1 MpcT 
(7) 



! ) (io : 



tot 
" \-iya 
9 erg s~ Mpc~ 

(8) 



6.6 



8.5 



+3.0 



1.1 



+0.6 
0.6 



-1.5(fix) 1.60 



1.1 



+0.9 
-0.8 



1.9 



+0.5 
0.4 



6.6 



+ 1.0 
0.8 



Note. — (1): Rcdshift. (2)-(4): Best-fit Schcchtcr parameters for 4>* and L L , respectively, a is fixed to —1.5. (5): Reduced 
X of the fitting. (6)-(7): Number densities and Lya luminosity densities calculated with the best-fit Schechtcr parameters down 
to the observed limit of Lya luminosity, i.e. log L-^ ya — 42.4 erg s~ . (8): Inferred total Lya luminosity densities integrated 
down to I/Lya = with the best-fit Schcchtcr parameters. 



with the best-fit Schechter function within the error bar 
that extends down to 0. This excess data point is solely 
made by one exceptional LAE that is an exten ded giant 
LAE, Himiko, reported in lOuchi et al.l ()2009al ). 

3.4. Evolution of Lya Luminosity Function 

Figure [7J compares our Lya LF of LAEs at z = 6.6 
with those at z = 3.1 and 5.7. Note that these low-z LFs 
are derived from la rge LAE samples of wide-field SXDS 
(|Ouchi et al.l 120081 ) with the same procedures (including 
cosmic variance errors) and similar data sets taken with 
the same instrument. In this sense, there are little sys- 
tematics among different rcdshift results caused by the 
sample selection and measurement technique. While the 
LFs do not change within error bars from z = 3.1 to 5.7 
(|Ouchi et al.ll2008l ). the LF decreases from z = 5.7 to 6.6 
beyond sizes of errors, i.e. uncertainties of statistics and 
cosmic variance. To quantify this decrease, we present in 
Figure |8] error ellipses of Schechter paramet ers of LFs at 
z = 6. 6, together with those at z = 5.7 from lOuchi et all 
(|2008f ). Figure [5] shows that error contours of z = 6.6 and 
5.7 differ at the > 90% significance level. Even if we con- 
sider the maximum effect of contamination (J J3.1| . there 
exist a significant evolution of Lya LF. Because a con- 
tamination correction pushes down our z = 6.6 LAE LF 
towards low number density in Figure [71 an inclusion of 
contamination correction even strengthens the evolution- 
ary tendency. The EWq limit of our z = 6.6 LAE sample 
is > 14A, which is different from and smaller than that 
of th e z = 5.7 LAE sa mple {EWq ^ 27A; lOuchi et all 
l200l by AEW ~ lOA. lOuchi et all (J2008I ) have inves- 
tigated a possible false evolution given by a choice of 
EWq limit for LAE samples at different redshifts. They 
use LAE samples at z = 3 — 6 whose EWq limits dif- 
fer by 30A(£W > 30 - 60A), and fit a number-^Wo 
distribution of the samples with a Gaussian function to 
obtain inferred Schechter parameters for all z = 3 — 6 
LAE samples with the same EWq limit of EWq = 0. 
lOuchi et al.l (|2008l ) find that the inferred </>* is different 
from the original cf>* only by < 10%, and claim that a 
choice of EWq limits provides a negligible impact on 
the results of Lya LF evolution. This is true if a dif- 
ference of EWq limit is much smaller than the best-fit 
Gaussian sigm a of number- EWq distribution (> 130A 
for z = 3 — 6; lOuchi et al.l 120081 ). Although we cannot 
carry out the similar investigation for our z = 6.6 sam- 
ple due to poorly constrained EWq values of our LAEs 
( ^6.1.11 for more details), it is very likely that the bias 
introduced by the different EWq limits (AfiWn ~ 10 A) 
is smaller than that found in lOuchi et "all (|2008f ). < 10% 



corresponding to < 0.04 dex in <j)* . This small bias of 
< 0.04 dex does not affect to the conclusion of LF de- 
crease in Figure EJ Moreover, the conclusion of LF de- 
crease is probably, again, even strengthened, because cor- 
rections for the incompleteness of the EWq limits raise 
ip* of the z = 5.7 sample {EWq > 27A) more than that 
of our z = 6.6 sample (EWq > 14A) in the case of the 
same number- EWq distribution at z = 5.7 and 6 . 6. We 
plot error contours obtained bv iKashikawa et all (J2006I ) 
in FigureElfor comparis on. Because the measurements of 
IKashikawa et al.l (J2006D do not include cosmic variance 
errors, their contours are relatively small. 

Although the errors of our measurements are not small, 
Figure [5] implies that a decrease in L* would be the dom- 
inant factor of the LF evolution from z = 5.7 to 6.6. We 
investigate whether pure number or pure luminosity evo- 
lutions are dominant. We fix either L* or (j>* of the z = 
5.7 LF, and carry out Schcchtcr function fitting to our 
z = 6.6 LF. Figure[7]plots the best-fit Schechter functions 
for the two scenarios of the pure number (dotted line) and 
pure luminosity (dashed line) evolutions from z = 5.7 to 
6.6. We obtain <j>*(z = 6.6) = 0.50* (2 = 5.7) for the pure 
number evolution, and L*(z = 6.6) = 0.7L*(z = 5.7) for 
the pure luminosity evolution. The fitting to the pure lu- 
minosity evolution is better than that to the pure number 
evolution, and the x 2 value of the pure luminosity evolu- 
tion is smaller than that of the pure number evolution by 
7.40. We confirm that, irrespective of a, the pure lumi- 
nosity evolution is more favorable than the pure number 
evolution based on our Schechter function fitting with 
the fixed best-fit L* or (j>* of the z = 5.7 LF in the cases 
of a = -1.0 and -2.0 (Table 5 of lOuchi et all 12001) 
that bracket a reasonable range of a. In the case of pure 
luminosity evolution, the decrease of luminosity is 30% 
from z = 5. 7 to 6.6. This decrease is smaller than that 
obtained by IKashikawa et al.l (|2006l ) who claim a lumi- 
nosity decrease by 0.75 magnitude. Moreover, this 30% 
decrease is too small to be distinguished by any previous 
studies that include large statistical and cosmic variance 
errors. We plot LF estimates from each of the five ~ 0.2 
deg 2 subficlds in Figure [7] with the same open symbols as 
in Figure [6] These open symbols illustrate that with the 
data of a single ~ 0.2 deg 2 field alone which is a typical 
survey size of previous studies it is difficult to distinguish 
whether or not z = 6.6 LFs show evolution (decrease) 
with respect to z = 5.7. We find the decrease of Lya LF 
between z = 5.7 an d 6.6 at the > 90% significa nce level. 
On the other hand, iMalhotra fe Rhoadsl (|2004t ) claim no 
evolution between z = 5.7 and 6.6, which appears to con- 
tradict with our finding. However, IMalhotra fc Rhoadsl 
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Fig. 9.— Sky distribution of the SXDS LAEs at z = 3.1 (left), 3.7 (center), and 5.7 (right) obtained by IQuchi et~aTI I2008T ). Red 
squares, magenta diamonds, and black circles present positions of narrow-band (bright, medium-bright, and faint) LAEs, respectively, in 
narrow-band magnitudes of (NB503 < 23.5,23.5 < AT.B503 < 24.0,24.0 < iVB503 < 25.3) (left panel), (NB570 < 23.5,23.5 < NB570 < 
24.0,24.0 < NB570 < 24.7) (center panel), and (7VB816 < 24.5,24.5 < AB816 < 25.0,25.0 < AB503 < 26.0) (right panel). The gray 
shades represent masked areas that are not used for sample selection. The scale on the map is marked in both degrees and the projected 
distance in comoving megaparsecs at each redshift. 

comes around the middle of previou s resul ts an d falls be- 
tween those of iMalhotra fc RhoadU (12004D and iHu et al.l 
(|2004) : the one used by IMalhotra k Rhoadsl (l2004f) is 
placed below ours, while the one used bv lHu et al.1 ( 20051 ) 
is above ours. Because we see exactly the same ten- 
dency in z = 6.6 L F measurements in Figure E] (i.e. 
IMalhotra fc Rhoadsl' s LF below ours, and lHu et all s LF 
above ours) , the 'no evolution' i s more clearly found in 
Figure 2 of IMalhotra k Rhoadsl (|2004f ) and Figure 3 of 
IHu et all (|2005l ). even without considering their large un- 
certainties of statistics and cosmic variance. 

In summary, we conclude that Lya LF decreases from 
z = 5.7 to 6.6, and that no evolution of Lya LF is ruled 
out at the > 90% significance level with errors including 
uncertainties of statistics and cosmic variance. The pure 
luminosity evolution is more preferable to the pure num- 
ber density evolution. If the pure luminosity evolution is 
assumed, L^ dims by 30% from z = 5.7 to 6.6. 
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Fig. 10.— Same as Figure [9] but for our 207 LAEs at 
z = 6.565 ± 0.054. Red squares, magenta diamonds, and black 
circles show positions of bright (7VB921 < 25.0), medium-bright 
(25.0 < NB921 < 25.5), and faint (25.5 < NB921 < 26.0) LAEs, 
respectively. The red square highlighted with a red open square 
indicates the giant LAE, Himiko, with a bright and e xtended Lya 
nebular at z = 6.595 reported by Ouchi et al. (20 09al) . 



(|2004D carefully state that a Lya flux of z = 6.5 LAEs is 
not attenuated by a facto r of 3, i.e. 300%, compared w ith 
z = 5.7 LAEs. Because IMalhotra k Rhoadsl (|2004f ) do 
not conclude evolutionary tendency in the change by a 
few 10% level , but by a few 100% l e vel. In this sense, the 
conclusion of IMalhotra k Rhoadsl (|2004t ) quantitatively 
agrees with ours. Another issue is the baseline LF at 
z = 5.7 used in these comparisons. The baseline z = 5.7 
LFs a re diffe rent betwe e n our study, IMalhotra k Rhoadsl 
(120041) . and IHu et al.l (120051) . We use the up-to-date 
5.7 LF of IQuchi et all (120081) sh own in Fi g ure T 



and IHu et al. 



while IMalhotra k Rhoadsl (J2004 
use their own LFs. Figure 19 of IQuchi et al. 
compares all of these three z = 5.7 LFs, and indicates 
that our baseline z = 5.7 LF of IQuchi et~aTI (|2008D 




4. SPATIAL DISTRIBUTION 

In this section, we discuss clustering of our z = 6.6 
LAEs. Since we need to know an evolutionary trend of 
clustering at z = 3 — 6 to distinguish cosmic reionization 
and galaxy evolution effects at z = 6.6 f £|6.1.2[) , we de- 
rive correlation functions of not only our 207 LAEs at 
Z = 6.6, but also those at lower redshifts. We use 356, 
101, and 401 LAEs at z = 3. 1, 3.7, and 5.7, respectively, 
given by I Ouchi et al.l ((2008). Since we have a moder- 
ately large number (> 200) of LAEs at z = 3.1, 5.7, and 
6.6 down to narrow-band magnitude limits of 25.3, 26.0, 
and 26.0, respectively, we make bright subsamples with 
narrow-band limiting magnitudes brighter by 0.5 magni- 
tude; 24.8 (z = 3.1), 25.5 (z = 5.7), and 25.5 {z = 6.6). 
Readers should refer IQuchi et al.l (|2008h for the details 
of the SXDS z = 3.1 - 5.7 LAE samples. 

4.1. Angular Correlation Function 

Figures IMTOl present sky distributions of the z = 
3.1 — 5.7 LAEs and our z = 6.6 LAEs, respectively. In or- 
der to quantitatively measure the inhomogeneity of spa- 
tial distribution, we derive the angular two-p oint correla- 
tion fu n ction (ACF) i n the s ame manner aslOuchi et al.1 
(I2003I) . IQuchi et all (l2004bl) a n d IQuchi et al.l (l2005bl). 
According to lLandv k Szalav I (|1993l ). the ACF is cal- 



Statistics of Lya Emitters at z ~ 7 



11 



Mpc h'j 
0.1 1.0 10 



Mpc hj 
0.1 1.0 10 



0.1 



Mpc h 7 ' 
1.0 10 



100 




10 100 1000 1 
9 (arcsec) 



10 100 1000 1 
9 (arcsec) 



10 100 1000 
9 (arcsec) 



Fig. 11. — Angular correlation function (ACF) and bias of the SXDS LAEs at z = 3.1 (left), 3.7 (center), and 5.7 (right). Top panels 
present ACFs, and bottom panels show bias as a function of angular distance. Black filled and open squares indicate the best estimates for 
all LAEs with and without the corrections of limited area (i.e. integral constraints) and number (eq. [8J. Gray filled and open squares are 
the same, but for the bright subsamples, NB503 < 24.7 in the left panels and NB816 < 25.5 in the right panels. Crosses in the top panels 
represent ACFs calculated with the Monte-Carlo simulation-based random sources that are artificial LAEs distributed and detected in the 
real images (see the text). In the top panels, solid lines are the best-fit power law funct ions for all LAEs, while the dotted lines denote 
ACFs of underlying dark matter predicted by the CDM model (Peacock & Dodds 1996). In the bottom panels, solid and dashed lines 
present the best estimates of bias and the assoc i ated la errors. In the top right panel, filled circles and dot-dashed line are ACF and the 
best-fit power law obtained by[Murayama et al. (2007). The scale on the top axis denotes the projected distance in comoving megaparsecs 
at each redshift. 

samples. We first create a pure random sample com- 
posed of 100,000 sources with the same geometrical con- 
straints as of the data s ample, and estim ate errors with 
the bootstrap technique (|Ling et al.lll986t ). Figures [TTIfT2l 
show the ACFs, cj bs(#), of LAEs from the observations 
at z = 3.1 — 5.7 and 6.6. We find significant clustering 
signals for our z = 6.6 LAEs as well as the z = 3.1 — 5.7 
LAEs. 

We then confirm that these clustering signals are not 
artifacts produced by the slight inhomogeneous quality 
over the images or occultation by foreground objects on 
the basis of our Monte Carlo simulations. We use mock 
catalo gs of LAEs obtained by simulations of lOuchi et al.l 
(|2008t ). which have number counts and color distribu- 
tion that agree with observational measurements. We 
generate 50,000 artificial LAEs based on the mock cat- 
alog, and distribute them randomly on the original 1 
deg 2 narrow- and broad-band images after adding Pois- 
son noise according to their brightness. Since most of 
LAEs are nearly point sources, we assume profiles of 
point-spread functions that are the same as the origi- 
nal images. Then, we detect these simulated LAEs and 
measure their brightness in the same manner as the real 
LAEs. We iterate this process 10 times and select LAEs 
with the same color criteria as the real LAEs. We thus 
obtain ~ 200, 000 simulation-based random sources at 
each redshift whose positions arc affected by the inhomo- 
geneity of LAE detectability and the occultation of fore- 
ground objects, and therefore slightly different from the 
pure random sample. We use these simulation-based ran- 
dom sources for our ACF calculation. Crosses in Figure 
[Tl]prcscnt estimates of ACFs with these random sources. 
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Fig. 12. — Same as Figurc[TT] but for our LAEs at z = 6.6. Gray 
filled and open squares are the estimates for the bright (NB921 < 
25.5) subsample. 

eulated by: 

w ohs (9) = [DD(0) - 2DR(9) + RR(8)} / 'RR{9) , (6) 

where DD(9), DR(9), and RR(9) are numbers of galaxy- 
galaxy, galaxy-random, and random-random pairs nor- 
malized by the total number of pairs in each of the three 
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TABLE 4 
Clustering Measurements 



2 


Sample 


C + l/N 


A u 

(arcsec^) 


/9 


ro 

(h- \ Mpc) 


„max 
1° 

Ciion M P C ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


3.1 


All (NB503 < 25.3) 


1.99 x 10 -2 


5.61 ±2.46 


0.80 (fix) 


J -' u -0.46 


I OQ+0.45 
1 - a! '-0.55 


3.1 


iV,B503 < 24.8 


1.66 x 10 -2 


3.13 ±3.12 


0.80 (fix) 


i 9 n+0.58 


1 44+0-67 
± - yi ^-1.37 


3.7 


All (NB570 < 24.7) 


5.59 x 10^2 


15.24 ±6.38 


0.80 (fix) 


9 yi+0.58 
z ''^-0.72 


o 9/L +0.69 


5.7 


All (JVB816 < 26.0) 


4.52 x 10 -2 


14.13 ±2.80 


0.80 (fix) 


o , 9 +0.33 
O,lz -0.36 


4 30+ 045 


5.7 


NBS16 < 25.5 


8.76 x 10~ 2 


26.37 ±6.59 


0.80 (fix) 


4 41+0-59 


6-08±g : |g 


6.6 


All (NB921 < 26.0) 


2.53 x 10 -2 


6.63 ±3.73 


0.80 (fix) 


9 O1+0.65 


3.60±i;i 


6.6 


NB921 < 25.5 


9.35 x 10~ 2 


23.25 ±22.31 


0.80 (fix) 


4 64+ 2 ' 10 


1 OQTU' ^O 

'• zo _6.01 



Note. — (1): Rcdshift, (2): LAE sample, (3): offset value for A w that corrects for limited survey area 
and number (cq. [8j, (4): best-fit amplitude of power-law function for the angular correlation function 
in arcsec". (5): the fixed /3 for the power-law fitting, (6)-(7): correlation lengths in /iTnn Mpc. The 
best-estimate and the maximally-contamination corrected ro in (6) and (7), respectively. To facilitate 
comparison with previous results, we express ro using hioOi the Hubble constant in units of 100 km s — 
Mpc -1 , instead of 70 km s _1 Mpc -1 . 



The ACFs estimated with these simulation-based ran- 
dom sources (crosses) and the pure random-distribution 
sources (squares) are consistent. Accordingly, we con- 
clude that the clustering signals are not artifacts given 
by the slight inhomogeneity of image qualities or the oc- 
cultation of foreground objects. 

To evaluate observational offsets included in w OD s(#) 
due to the limited area and object number, we assume 
the real ACF, uj(9), is approximated by the power law, 

u(6) = A u 9-P. (7) 

Then, the offset from the observ ed ACF, uj i, s (9), is give n 
by the integral constraint, C, (jGroth k Peebles! Il977[ ). 
and the number of objects in the sample, N . 

uj(6)=uj obs (6) + C+±- (8) 



N 



C = 



Y>RR(9)A U 9~ 



(9) 



Y,RR{9) 

The term, 1/N, in cq. (JSj) corrects for the difference 
between the number o f object pairs, N(N — l)/2, and its 
approximation, N 2 /2 ((Peebles 1119801 ). Note that most of 
the previous clusterin g studies for high-z ga l axies neglect 
this 1/N term (e.g. JRoche fc Ealesl 119991 : Ibaddi et al.1 
|2000( ) probably because of their large samples. However, 
it should be applied for samples with a small number of 
objects such as LAE samples to obtain more accurate 
ACF at a large scale. The ACF corrected with eq. (JSJ 
are also presented in Figures [TT1fT2l . 

We fit the power law (eq. 0) over 10" < 6 < 1000" 
with the corrections. The lower limit of the fitting range, 
9 = 10", is placed, because 1-halo t erm of high-z galax- 
ies is dominant at this small scale (lOuchi et alj|2005bt 
iLee et al.|[200l Iffildebrandt et al.ll2007ll2009t) . The up- 
per limit of the range is nearly the limit of our ACF mea- 
surements with significant signals. We apply the upper 
limit of 400" for the z = 5.7 LAEs, since a large-scale 
clustering (> 20 Mpc ) with proto-cluste rs is reported 
in this survey volume (jOuchi et al.ll2005al ). Because we 
obtain no meaningful constraints on /3, we fix (3 with 
the fiducial value of 0.8 following t he previous cluster- 
ing analyses (e.g. lOuchi et alJl2003t iGawiser et al J 120071 
iKovac et al.ll2007l ). We summarize the best-fit values in 
Table |H 



Foreground contamination to a galaxy sample dilutes 
the apparent clustering amplitude of galaxies. If the frac- 
tion of contaminants is f c , the apparent A u value can be 
reduced by a factor of up to (1 — f c ) 2 . The correlation 
amplitude with the contamination correction, A™ ax , is 
given by 

*r - (T^f do) 

This is the maximum reduction of the correlation am- 
plitude that occurs when the contaminants are not at 
all clustered. In reality, contaminants, if any, are the 
sum of foreground galaxies mostly at some specific red- 
shifts, and thus would be clustered to some extent on 
the sky. We use the maximum values of / c , and obtain 
maximal A Illax values inferred from the LAE data. The 
maximum value s of f c are (0.13,0 .14,0.25) for redshifts 
of (3.1,3.7,5.7) (lOuchi et al J 120081) and 0.33 for z = 6.6 
f £|3.1[) . In the next section, we will place conservative up- 
per limits on correlation lengths and bias measurements 
with the maximal values of A™ ax . 

4.2. Correlation Length and Bias 

The spatial correlation function of galaxies is well- 
approximated by a power law as 

£ = (r/ro)"T (11) 

only with a sub tle departure from the real spatial corre- 
lation function (jZehavi et al.l 12004 ). where r is the spa- 
tial separation between two objects, ro is the correlation 
length, and 7 is the slope of the power law. The corre- 
lation length, ro, is related to the correlation amplitude, 
A u , with the integral equation called Limber equation 
([Peebles 1119871 lEfstathiou et al.|[l99lT) . 



A 



= kr a 



F(z)Dl-'(z)N(z) 2 g(z)dz 



N(z)dz 



(12) 

where F(z) 3 describes the rcdshift dependence of £(r), 
Dg(z) is the angular diameter distance, N(z) is the red- 



lxi Assuming that the clustering pattern is fixed in comoving 
coordinates in the redshift range of our sample, we take the func- 
tional form, F(z) = (1 + z)/(l + 2 c ) - ' 3+e ', where z c is the central 
rcdshift of the sample and e = —1.2. The effect of the change in e 
over < e < —3 on ro is, however, very small. 
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TABLE 5 
Bias and Hosting Dark Halos 



z log(L LiE 
(erg s -1 
(1) (2) 



(Mpc" 3 ) 

(3) 



(1) 



(5) 



M h 
(6) 



jymin 

(7) 



1 / max 

(M e ) 
(8) 



3.1 


42.3 


8-9ir. 9 2 


3.1 


42.1 


l-5± U n1 


3.7 


42.6 


9 q+0.6 


5.7 


42.6 


3.5±?| 


5.7 


42.4 


6.8ii;? 


6.6 


42.6 


i a+0-5 


6.6 


42.4 


4 1+0.9 

^-0.8 



x 10- 
x W 
x 10- 
x 10- 
x 10- 
x 10- 
x 10- 



1.5 ±0.7 
1.5 ±0.3 
2.8 ±0.5 
6.1 ±0.7 

5.5 ±0.4 
6.0 ±2.2 

3.6 ±0.7 



1.7 ±0.8 

1.8 ±0.3 
3.2 ±0.5 

8.2 ±0.9 

7.3 ±0.5 

8.9 ±3.3 

5.4 ±1.0 
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Note. — (1): Rcdshift, (2): approximated limiting Lya luminosity estimated from narrow-band magnitudes, (3): 
number density of the observed LAEs that is derived with the best-fit Schcchtcr function ( Table [3l, (4): best estimate of 
bias, (5): upper limit of bias with the maximal contamination correction, (6): average hosting dark halo mass inferred 
from the best estimate of bias, (7): minimum hosting dark halo mass, and (8): upper limit of the average hosting dark 
halo mass corresponding to the upper limit of bias. 



shift distribution of objects, g(z) = Hq/c[(1 4- z) 2 (l + 
Q-mZ + £|a.((1 + z)~ 2 — l)) 1 / 2 ], and fc is a numerical con- 
stant, k = V7tT[(7 - 1 )/ 2 ]/ r (7/ 2 )- Thc slo P e > P, of thc 
angular correlation function is related to 7 by 

7 = /? + !• (13) 

We adopt the redshift distribution of the LAEs calcu- 
lated with the simulation output sources (i.e. referred to 
as simulation-based random sources in our Monte Carlo 
simulations; §4.1j) . Thc correlation lengths thus obtained 
are summarized in Table SJ These correlation lengths are 
expressed with /iioo, the Hubble constant in units of 100 



km 



instead of 70 km s , to facilitate comparison 



with previous results. In Table |4j we present tq and 
its upper limit, r™ ax , that are the correlation lengths 



obtained from A u , and its maximal value, A" 



respec- 



tively. The correlation lengths are typically 1 — 2 h 100 
Mpc at z ~ 3 and 2 — 6 /i7/ q Mpc at z = 4 — 6. Our 
estimat es of rp are consiste nt wit h the previous mea sure- 
ments. lOuchi et all (|2003l ) and iKovac et all (|2007t ) ob- 
tain (r , rff**) of (3.5±0.3, 6.2±0.5) h^ Mpc at z = 4.9 
and (3 .20±0.42,4.61±0.60) ftJ Mpc at z = 4 . 5, resp ec- 
tively. iGawiser et all (|2007f) and lGuaita et al.l (|2010[ ) es- 
timate r to be 2.5t°;7 ^ioo M P C at z = 3.1 and 3.4 ±0.6 
/i7/oo Mpc at z = 2.1, respectively, where r is converted 
to values in units of feTnn Mpc. At z = 3.1, the value of 
IGawiser et al.l (|2007t ) is slightly higher than that of our 
z = 3.1 LAEs, but these are still consistent within la 
errors. The correlation length of our z = 6.6 LAEs is 



'•() 



2 — 5 h wo Mpc, while that for 



„max 




7 h 



100 



Mpc, as presented in Table |U 

Note that these correlation leng ths are estim ated by 
the Limber equation of cq. (|12p . ISimonl (|2007[ ) claims 
that the Limber equation is accurate for small galaxy 
separation but breaks down beyond a certain separation 
that mainly depends on a width of galaxy distribution 
and distance. In this sense, eq. (|12p is regard as Lim- 
ber approxi mation t hat is valid only in a small angular 
separation. ISimonl (J2007I ) presents the exact equation 
connecting the spatial and an gular co rrelation functions, 
which is he reafter referred as iSimonl 's equation. On the 
other hand, iSobral et al.l (|2010f ) have demonstrated that 
there is only a small, ~ 4%-level, difference bet ween th e 
correlation lengths calculated with Limber's and ISimonl 's 



equations based on their NB-selected sample of Ha emit- 
ters at z = 0.8. Although their survey area is similar to 
that of our study, it is not clear whether the Limber ap- 
proximation gives only a negligible impact on our results 
at the different redshift of z = 6.6. We hav e compared 
ACFs calculated with Limber's and ISimon|'s equations 
in the same manner as shown in Figure 4 of ISobral et al.l 
(|2010l) . We have found that there are no apparent differ- 
ences between two results up to 9 ~ 200" . We have used 
our data of 10" — 1000" for our power-law fitting that is 
mostly not affected by the bias of Lim ber app roximation. 
Moreover, the ACFs of Limber's and ISimonl 's equations 
differ only by a factor of ~ 2 even at 9 near 1000" , the 
largest angular separation of our fitting where the er- 
ror of our data point is comparably large. Thus, the r$ 
values change negligibly small for our galaxies at higher 
redshifts, compared with their relatively large statistical 
errors. 

We define galaxy-dark matter bias in the framework of 
ACDM model by 



b g {9f=uj{9)/uj dm {9), 



(14) 



where uiam(O) is the ACF pred icted by the non-linear 
model ([Peacock fc Dod ds 1996) at the cosmic volume 
defined with N(z) for each sample. Note that Wdm(^) 
is defined in the same survey volume as our LAE sam- 
ple, and that this bias value is equivalent to bias given by 
the original definition, b^ = £(r) / ^ m (r) , where £dm(0 is 
the spatial correlation function of dark matter. Top pan- 
els of Figures [TTiri"2l plot Wd m (fl), and the bottom panels 
of Figures ITT1fr2l show b s (9). We obtain an average bias, 
6 g , that is a value of b g (9) averaged over the same angu- 
lar range as the one for the A^ fitting ( £14. ip with error 
weights. The upper limits of bias values, &™ ax , are also 
calculated with A" lax for thc maximal contamination cor- 
rection case of eq. (fTU| . The estimates of b g and 6™ ax are 
presented in Table [5] For our z = 6.6 LAEs, we obtain 
b e = 3 - 6 and ^ lax = 5-9. 

4.3. Hosting Dark Halo 

We estimate dark halo masses of LAEs with bias val- 
ues obtained in Section 14.21 in the framework of ACDM 
model. Because, as discussed in ^6.1. 2[ we find that clus- 
tering of z = 6.6 LAEs is negligibly affected by cosmic 
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Fig. 13. — Bias of LAEs at z = 2 — 7. Red stars present our LAEs 
at z = 6.6 with the bright (top) and all (bottom) samples. Cyan, 
blue, magenta stars plot bias of SXDS LAEs at z = 3.1, 3 . 7, an d 
5.7, respectively, estimated with the samples ofRjuchi et al. (120081 '! . 
Small offsets are given to these points in the direction of redshift to 
avoid overlapping symbols. Right and Left cyan stars denote the 
bright and all samples at z = 3.1, respectively. For z = 5.7 LAEs, 
we show the bias of the bright LAE sample at right, all LAE sample 
in the m iddle, and the similar sample but estimated bv lOuchi et al.l 
(2005a) (O05) with the count-in-cell technique in a large scale of 
20-Mpc radius at l eft. Black stars represent estimates g iven by 
IGuaita et a"L1 IMrTi (G u lO; z = 2.1), IGawiser et al.l 1200711 (Ga07; 
z = 3 .1), [Kovac ct aL| J2007l1 (K07; z = 4.5), and IQuchi et~aTI 
l2003h (O03; z = 4.9 ). Because the fields of the SXDS LAEs at z = 
5.7 and lOuchi et a l. ( 2003) have a strong large-scale clustering with 
proto-cluster candidates made of LAE overdensity, we use open- 
star symbols to indicate the existence of these structures that may 
boost the bias measurements. Solid lines denote bias as a function 
of redshift for dark halos with a mass of 10 s , 10 9 , 10 10 , 10 11 , and 
1O 12 M predicted by the ISheth fc Tormenl l Tl99Sf) model in the 
case of one-to-one correspondence between galaxies and dark halos. 
Gray area denotes the dark halo mass range of 10 10 — 10 12 AfQ 
where LAEs at z = 2 — 7 falls. Dotted lines represent evolutionary 
tracks of bias under the galaxy-conserving model (eq. I24[ i. This 
plot assumes as = 0.8. We apply a correction to all o f the previous 
estima tes that assume erg = 0.9, except for the one of lGuaita et al.l 
(IMIiih . 

rcionization, compared with the size of errors, bias of 
Z = 6.6 LAEs mostly depends on properties of hosting 
dark halos. 

Generally, the average dark halo mass, Mh, and the 
mean galaxy bias, b g , are given by 



M h = 



ggn M ft$(M) N g (M) n(M) dM 

$M^ f^(M) N e (M) n(M) dM 
J~„ in b(M) ft$(M) N g (M) n(M) dM 
fur* CM N e (M) n(M) dM 



anfll5) 



(16) 



where n(M) and b{M) arc number density and bias of 
dark halos, respectively, with a mass of M, and N g (M) 
is a galaxy occupation function at the halo mass, M. 



to 



/duty (-^0 i s a fraction of dark halos hosting LAEs 
those hosting any galaxies, which we refer to as duty cycle 
of LAEs. M h nm is the minimum dark halo mass that can 
host a g alaxy. To obtain n(M) a nd b(M), we apply the 
mode l of ISheth fe Tormenl (j!999f ) (see also iMo fe White! 
2002). 



n{M)dM = A{1 



1 2 po dv' ( v 12 



.I2q 



■ cxp 



7T M dM """ V 2 

(17) 
where v' = y/av, a = 0.707, A ~ 0.322 and q = 0.3. p 
is the current mean density of the Universe. Here, v is a 
function of the growth factor, D(z), and the rms of the 
density fluctuations on mass scale M, <r(M), as defined 
by 

VS D(z)*(M) ' (18) 

where S c = 1.69 represents the critical ampli- 
tude of the perturbation for collapse. We calcu- 
late the_growthi_^£to£ ! __^z) ! _at redshift, z, follow- 
ing iCarroll. Press, fc Turner! (|1992[ ). Again the power 
spectrum of the densi ty fluctuations is calculated with 
the tra nsfer function of lBardeen. Bond. Kaiser, fc Szalavl 
(|1986t) . The bias of dark halos is estimated by 



1 
To 



v a + W 2{x-c) 



:c + &(l-c)(l-c/2) 



(19) 

where b = 0.5, c = 0.6, and the other parameters (z/, S c , 
and a) are the same as those in eqs. (fl"7| and (p~8|) . 

Note that our clustering measurements do not have ac- 
curacies high enough to determine the galaxy occupation 
function which is often described with three free param- 
eters (see e.g. JHamana et al.l 12004 IQuchi et al.l [2005b; 
iLee et al.ll2006l : iKovac et al.ll2007l ) Different choices of 
galaxy occupation functions change average halo masses 
only by a factor of a few, but not ove r an order of mag- 
nitude for LAEs (JHamana et al.lr2004f ). unless extremely 
unphysical parameter sets are taken. Because we aim 
to accomplish an accuracy of an order of magnitude 
for halo properties inferred from clustering, we assume 
N g (M) = 1, a one-to-one correspondence between galax- 
ies and dark halos. We also assume that duty cycle of 
LAEs is mass independent for simplicity (see more dis- 
cussions about duty cycle in §4.4|) . In this case, Mh and 
6 g are irrespective of /^ty ■ We determine M™ m with 
our observational estimates of LAE bias, b g , via eq. (|l~6"|) . 
and calculate the average dark halo mass, Mh, from the 
M^ in values with eq. ([IS). Table [5J presents M h and 
M™ of our LAEs, together with the maximum of av- 
erage halo mass, M" lax , estimated with the maximally 
contamination corrected bias of 6 max . 

Figure [13] plots our b g values as a function of red- 
shift. This figure also compares b g of LAEs obtained 
in the previous studies. Note that in our study we as- 
sume us = 0-8 that is consistent with the latest W MAP 
results (jKomatsu et all 120091 : lLarson et ail 120101) . Be- 
cause all of the previous studies, except IGuaita et al.1 
(|2010l ). assume cr% = 0.9, we multiply the bias val- 
ues of these previous results by 1.1 to correct for the 
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Fig. 14. — Bias as a function of number density at z ~ 3 (left), 2 ~ 4 (center), and jz ~ 6 — 7 (right). Star marks represent LAEs, 
while the other symbols denote LBGs (including dropout galaxies). Specifically, the red stars indicate our LAEs at z = 6.6. The open 
star marks in the rig ht panel are z = 5.7 SXDS LAEs whose bias values might be boosted by the large-scale c lu sterin g with proto-clusters 
llOuchi ct al. 20053). All of the star marks are our estimates from samples of this study and IQuchi et alT s 12003) studi es, except the 
small star mark in the left panel that is given by Gawiscr ct al. (2007). Because our estimates of LAEs overlap with that of lGawiser et all 
l)2007T ). we sligh tly shift their point in the number densi ty direction for clarification. We plo t the bias and numbe r densities of LBGs with 
squ ares (z ~ 4; [Ouchi et^d.| l2004bl) , diamonds (z ~ 4; | Lee et at) 120061 ), pent agons (z ~ 6; IQver zier ct al.N2006l) , hexagons (z ~ 3.3 and 
3.8; IHildebrandt et al.1120091) , and circles (z ~ 4 IQuchi et ~ . 2005b and z ~ 6; IQuchi fc SxdsT cam 1120051) )■ Bl ack solid lines indicate the 
relation of bias and number density for the case of one-to-one correspondence between galaxies and dark halos (Shcth & Tormcn 1999) at 
z = 3, 4, and 6, in the left, center, and right panels, respectively. Black dashed lines present the same relation but their number densities 
are multiplied by 1/10 and 10 from bottom to top. Red solid and dashed lines in the right panel are the same, but for z = 6.6, the redshift of 
our LAE samples. This plot assumes erg = 0.8. Because all of the previous bias estimates here, are obtained with a cosmological parameter 
of erg = 0.9, we apply the correction to the previous estimates for the different erg values. 

difference of <j s . This correction factor of 1.1, irre- 
spective of redshift and selection function, is estimated 
from y/0Jdm(crx = .9)/Lj^ m (^H = 0-8) with the model of 
iPeacock fc Doddsl (|1996l) . where Wd m (o'8 = 0.9) and 
Wdm(o'8 = 0.8) are ACFs predicted by the non-linear 
model under a$ = 0.9 and 0.8, respectively. Figure [13] 
shows that our all and bright subsample measurements 
{z = 3.1 and 6.6) agree within la errors at each redshift. 
This indicates that a clear luminosity segregation cannot 
be found beyond the size of our relatively large error bars 
for samples whose luminosity differs only by a factor of 
~ 1.6 (see Tables [4j5]). The LAE samples of the previ- 
ous studies and ours have similar Lya luminosity limits 
(Lhya — 2 x 10 42 erg s _1 ) only with small differences 
by a factor of 2, and the effect of luminosity dependence 
is probably smaller than the given relatively large error 
bars. Thus, we consider that the luminosity dependence 
between the LAE samples is negligibly small in Figure 
[T3l (see also the arguments for the relation between bias 
and number density described below). Figure 1131 shows 
that our measurements at z = 3.1 arc consistent with the 



one at the same redshift given by iGawiser et al.l ()2007| ) . 
The bias of LAEs appear to increase from z = 2 — 3 
to 6.6. Because there are strong large-scale clustering 
with proto-clusters or proto-cluster ca ndidates reported 
in th e data of z = 4.9 in the S DF (ISh imasak u et al.l 
I2003T ) and z = 5.7 in the SXDS (jOuchi et all l2005al) . 
Figure [TBI plots open st ar marks for the b ias estimates 
at Z = 4.9 in the SDF dOtichi et al ] [20031) and z = 5.7 
in the SXDS (jOuchi et al. 2005a| ) to indicate that these 
bias estimates might be boosted by these obvious cosmic 
structures. Even if we omit the points of these open star 
marks in Figure 1131 we confirm the trend of increasing 
bia s toward high-z. T o under stand this increase, we plot 
the ISheth fc Torment s (|1999l ) model of a constant dark 
halo mass under the assumption of one-to-one correspon- 
dence between galaxies and dark halos. The dark halo 
masses of our z = 6.6 LAEs are about 10 10 — 10 n Mo 
in Figure [131 which are the same as those listed in Ta- 
ble [5] At z = 2 — 7, the bias values fall in the range of 
~ 10 10 — 1O 12 M for dark halo masses. Thus, the aver- 
age dark halo mass of LAEs is roughly ~ lO n±1 M0 at 
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z = 2 — 7. The dark halo masses of LAEs show no sig- 
nificant evolution at z = 2 — 7 beyond the mass-estimate 
scattering by an order of magnitude. It would indicate 
that LAEs, observed so far, are galaxies at the evolu- 
tionary stage for all or some type of galaxies whose dark 
halos have reached a mass of ~ 1Q 11±1 Mq- 

Figure [TJ] plots bias and number density of LAEs, and 
compares those with Lyman break galaxies (LBGs) in- 
cluding i-dropouts. The solid lines present the relation of 
bias and number density for all dark halos, which corre- 
spond to the case of one-to-one correspondence between 
galaxies and dark halos. The dashed lines are the same 
as the solid lines but with their number densities mul- 
tiplied by 10 and 1/10. The bias values of LAEs are 
generally smaller than those of LBGs at redshifts over 
3 — 7. One exception is LAEs at z = 5.7 that have a 
bias value as large as those of z ~ 6 LBGs. However, 
bias measurements of these LAEs could be boosted by 
strong large-scale clustering with proto-clusters as dis- 
cussed above. On the other hand, LBGs have a roughly 
one-to-one correspondence between dark halos and galax- 
ies within a scatter of one-order of magnitude in number 
density. Comparing with LBGs, we find that LAEs have 
significantly less number density and/or bias. This weak 
bias at the small number density for LAEs indicates that 
LAEs reside in less massive dark halos on average, and 
that not all but some of less massive dark halos can host 
LAEs. Because it is unlikely that some dark halos always 
have LAEs and the other dark halos with the same mass 
never host LAEs, dark halos become the ones hosting 
LAEs by duty cycle, stochastic processes made by star- 
formation history and changes of interstellar m edium ge- 
ometry and dynamics (Nagaminc et al. 2008). 

4.4. Duty Cycle 

In Figure HH the simplest estimate of duty cycle can be 
obtained from the ratio of the number density (ordinate 
of Figure Q3]) of all halos (solid line) to that of LAEs 
(star marks). However, due to a mild gradient of the solid 
curves, especially in the small bias range, even very small 
uncertainties of bias estimates of LAEs give very different 
duty cycle values. By this reason, a determination of 
duty cycle requires a precision measurement of bias of 
LAEs. Nevertheless, in the following subsection, we try 
to constrain the duty cycle of LAEs with the given bias 
measurements and errors of the present study. 

The duty cycle of LAEs can be constrained with a 
galaxy number density, n g , 



M° 



f^(M) N S (M) n{M) dM. (20) 



Following the assumptions, N E (M) = 1 and constant 



/du1v> ma de ii 



'duty 

expressed as /^ae _ «g//~ 



the equation of cq. (|20[) can be 



i(M) dM. We esti- 
mate /duty with M™ m values and the observed num- 
ber densities of LAEs (Table |5J to be f^f = (0.008 ± 
0.03, 0.04 ± 0.1, 1 ± 1, 0.009 ± 0.04) for our all LAEs at 
z = (3.1,3.7,5.7,6.6). Note that errors of f^S are raised 
by the combination of M" lm and n g estimates. The un- 
certainties of our ,/^y estimates are quite large as dis- 
cussed above. Especially, the /^y value of z = 5.7 pro- 
vides no meaningful constraints, because its error of n g 



is larger than the others due to a relatively large cosmic 
variance observed in our survey volume (see the caption 
of Figure HJ. Because of these large uncertainties, we 
cannot constrain the evolution of duty cycle. However, 
these /j^y estimates indicate that / dln ^y is a few 0.1 to 
a few percent, roughly ~ 1%. 

By definition, the duty cycle of LAEs, f^^B, is a frac- 
tion of galaxies that are LAEs. It can be expressed by 
the product, / d L u ^ = /|F^ x fg° t whcre /d SF y is a duty 

cycle of star-forming activities in all dark halos, and 
/duty ^ s * ne one for Lya emitting galaxies among star- 
forming galaxies. I Ouchi et al.l (|2004b[ ) estimate /fu ty of 
LBGs, UV-bright star-forming galaxies, to be 0.1 — 1 
from their clustering analyses, and a fainter galaxy 



have a smaller f^* t of 0.2 



at i 



25.5. Similarly, 



iLee et all (J2009D es timate /gf tv to be 0.15-0.6 for their 
LBGs with z' < 27 (|Lee et all 1 20061 ). Because our LAEs 
have a UV continuum magnitude of ~ 27 mag o n aver - 
age (|Orio et"aLll2010allbh . results of lOuchi et all (|2004b[ ) 
and lLee et all (J2009D suggest that /|,f ty is < 0.2 for LBGs 
whose UV luminosity is comparable to our LAEs. On the 
other hand. lStark et all (|2010[ ) obtain f^" y ~ 0.5 at - 27 
mag, which is a fraction of Lya emitting dropouts with 
EWq > 50A to all LBGs at z ~ 3 - 6. Combining these 
previous results, we find / d L u ^ = f^ x ffi" < p.i 
which is consistent with our estimates. lHaves et al.l 
(|2010f ) claim f d ^Z = 0.05 from the comparison of Lya 
and Ho: emitting galaxies based on Lya and Ha luminos- 

•SF ( 
duty 



ity functions. If we naively assume that /J 5 ,^.,, of Ha emit 



ting galaxies is the same as that of LBGs, i.e. f^ ty < 0.2, 
we obtain /^e < 0.01 (= 0.2 x 0.05), which is roughly 
consistent with our result. 

5. LYa LINE PROFILE 

We investigate Lya line profiles of our z = 6.6 LAEs, 
exploiting our high quality medium-high resolution spec- 
tra of Keck/DEIMOS. In this analysis, we add the 3 du- 
plicate spectra to the 16 spectra f £)2.3[) , and use a total 
of 19 DEIMOS spectra for z = 6.6 LAEs. In order to 
study the evolution of Lya line profiles between z = 5.7 
and 6.6, we also perform analysis of 11 DEIMOS spectra 
of z = 5.7 LAEs. Ten o ut of the 11 DEIMOS sp ectra 
are taken in the SDF by iShimasaku et al.l (|2006f ) with 
exactly the same instrumental configuration as that of 
our z = 6.6 spectra, and the remaining one DEIMOS 
spectrum is obtain ed for the SXDS z = 5.7 LAE sample 
(jOuchiet all 120081) as mask fillers of our DEIMOS ob- 
servations on 2008 October 3 f £|2.3p . The top panels of 
Figure [TBI present all of the DEIMOS spectra at z = 6.6 
and 5.7. 

5.1. Composite Spectra and Lya Velocity Widths 

We make composite spectra of our 19 spectra of z = 6.6 
LAEs, and the 11 spectra of z — 5.7 LAEs to inves- 
tigate the average Lya line profiles. The line centers 
of the spectra are aligned to the wavelengths of 9190A 
and 8157A that are average ones of the z = 6.6 and 5.7 
samples, respectively. The line fluxes arc normalized to 

14 These numbe rs are given in the sixth column in Table 6 of 
lOuchi et all <2004bl) . 
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Fig. 15. — Lya lines of LAEs. Top left: Red line represents the 
composite spectrum of our z = 6.6 LAEs. Black lines denote our 
19 DEIMOS spectra whose Lya lines arc aligned to the average 
redshift and scaled to the flux of the composite spectrum. Bottom 
left: The composite spectrum of our z = 6.6 LAEs is shown with 
a thick red line. One sigma error of the composite spectra are 
shown with thin red lines. Black curve presents the instrumental 
resolution of our DEIMOS spectra whose shape is approximated 
with a Gaussian profile. Top right: Same as the top left panel, but 
for the 2 = 5.7 LAEs. Blue line is the composite spectrum of the 
z = 5.7 LAEs. Bottom right: Same as the bottom left panel, but 
for the z = 5.7 LAEs. Blue lines are used for z = 5.7 LAEs. 

average Lya fluxes, since we calculate rejected-mean and 
median spectra for our stacking. In this process, we re- 
move lowest-quality five and two spectra from the z = 6.6 
and 5.7 LAE spectra, respectively, so as to avoid possi- 
ble systematic bias given by uncertainties of line center 
determination and/or flux normalization. We compute 
rejected-mean spectra weig hted with S/N rati os of Lya 
lines, in the same manner as lOuchi et al.1 (j2008f ). At each 
wavelength bin, the rejection is made for data points that 
fall beyond 68 percentile among all data points, which 
corresponds to the rejection of two (one) data points 
(point) above and below the 68 percentile for the LAEs at 
z = 6.6 (5.7). By this rejection algorithm, our compos- 
ite spectra are secure against systematic residuals of sky 
subtraction. We compare these composite spectra with 
the ones computed by simple median statistics, and check 
whether these rejection procedures introduce a system- 
atic bias. We find that the two composite spectra agree 
very well and that the rejected-mean composite spectra 
do not show any sign of bias in the overall shapes and ve- 
locity width of Lya lines. Because improvements of S/N 
are generally better in the rejected-mean stacking than 
in the simple median stacking, we use the rejected-mean 
spectra for our composite spectra in the following discus- 
sion. Figure [TBI plots our composite spectra of z = 6.6 
and 5.7 LAEs. The quality of the spectrum of z = 5.7 
LAE is not as good as that of z = 6.6 LAE. This is be- 
cause i) the total number of available z = 5.7 spectra 
(11) is small and ii) the average S/N ratio of individ- 
ual z = 5.7 spectra is rel atively low due to their short 
(~ 2 hr) integration time (IShimasaku et al1l2006lh The 



Fig. 16.— FWHM velocity width of our z = 6.6 LAEs, together 
with the 2 = 5.7 LAEs. Red square and blue circle denote velocity 
widths of the composite spectra. Black squares and open circles are 
velocity widths measured for individual LAEs at z = 6.6 and 5.7, 
respectively. Dotted and dashed lines represent the best-fit linear 
functions to the individual LAEs at z = 6.6 with and without mea- 
surement errors in the fitting, respectively. These FWHM velocity 
widths are the ones given by the Gaussian profile fitting, which 
provides robust measurements for relatively poor quality data of 
individual spectra. 



instrumental spectral resolution of the DEIMOS data are 
presented in the bottom panels of Figure [TBJ which illus- 
trate that our composite spectra are well resolved beyond 
the instrumental resolution. Our composite spectra show 
a very clear asymmetric profile because of the high S/N 
ratio. 

To quantify the Lya line profiles, we evaluate an 
FWHM velocity width, AVfwhm- We fit a Gaussian 
profile to the individual spectra as well as the com- 
posite spectra, and obtain AVfwhm- We correct for 
the instrumental broadening of line profile, and esti- 
mate intrinsic AVfwhm by AVfwhm = v^'obs = w mst> 
where i> b s and t>i nst arc FWHM velocity widths for the 
measured Lya lines and the instrumental resolution, re- 
spectively. We estimate errors of AVfwhm including 
uncertainties raised by a choice of profile fitting range 
and Poisson statistics based on Monte-Carlo simulations. 
Figure [16] and Table [2] present the AVfwhm measured 
from the spectra and total Lya luminosity estimated 
from the spectroscopic redshifts and the imaging data. 
Figure [TBI indicates that the AVfwhm values of z = 6.6 
LAEs are distributed in the range of 100-400 km s _1 at 
logi(Lya) ~ 42.6 - 43.6. The data points of z = 5.7 
LAEs are similarly distributed, but with a larger scatter 
given by the larger errors than those of z = 6.6 LAEs. 
Our composite spectra have the FWHM velocity widths 
of 270 ± 16 and 265 ± 37 km s" 1 at Z = 6.6 and 5.7, 
respectively. We also evaluate additional possible errors 
of AVfwhm introduced by uncertainties of line-center 
alignment in the process of spectrum stacking. We find 
that the typical line-center determinations are as good 
as 0.32A and 0.90A for z = 6.6 and 5.7 LAE spectra, 
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Fig. 17. — Evolution of LAE spectra. Red and black thick lines 
represent composite spectra of LAEs at z = 6.6 and 5.7, respec- 
tively. For comparison of line shapes, the composite spectrum of 
the z = 5.7 LAEs is redshifted to z = 6.6, and scaled by an arbi- 
trary factor. One sigma errors of the composite spectra are shown 
with thin red lines for z = 6.6 LAEs and a gray shade for z = 5.7 
LAEs. 

respectively, and that the additional errors of line-center 
alignment contributing to AVfwhm arc negligibly small, 
only a 0.1-1% level. We compare AVfwhm of the z = 6.6 
and 5.7 composite spectra in Figure [16] We find that 
AVfwhm does not evolve from z = 6.6 to 5.7 beyond the 
la level. Because the Lya luminosities of our compos- 
ite spectra are almost the same, logL(Lya) ~ 43.0, this 
comparison includes little bias of luminosity dependence. 

Note that these AVfwhm are estimated by Gaussian 
profile fitting. The real Lya profile is not Gaussian, but 
asymmetric, as clearly found in our composite spectra. 
Although one cannot obtain AVfwhm without fitting of 
a profile, such as Gaussian, for our relatively poor S/N 
data of individual spectra, AVfwhm can be directly mea- 
sured at least for the high quality composite spectra. The 
direct measurements of FWHM velocity widths corrected 
for the instrumental broadening are 251 and 260 km s _1 
for the z = 6.6 and 5.7 composite spectra, respectively. 
Table [6] summarizes the properties of our composite spec- 
tra, and indicates that the FWHM velocity widths of 
the direct measurements are comparable to those of the 
Gaussian fitting within the errors. 

Figure [T7] compares the composite spectra of z = 6.6 
and 5.7 LAEs. The composite spectrum of z = 6.6 LAEs 
is similar to the one of z = 5.7 LAEs within the errors. 
We find no large evolution beyond the la errors of our 
measurements. In fact, as shown in Table [HI the direct 
(Gaussian) FWHM velocity widths obtained above are 
251 (270) ±16 km s" 1 at z = 6.6 and 260 (265) ±37 km 
s _1 at z = 5.7, and there are no differences between the 
z = 6.6 and 5.7 composite spectra beyond the la error of 
~ 40 km s _1 , which is dominated by the er ror of z = 5.7 
composite spectrum. iHu fc Cowid (J2006I ) also find no 
evolution of line profiles, although there are no quantita- 
tive comparisons in their study. Our composite spectrum 



of z = 5.7 has an S/N ratio lower than that of z = 6.6 
by a factor of 2 in the FWHM measurements, which do 
not allow us to investigate a slight profile change from 
z = 6.6 to 5.7 that may exist. If one compares the best- 
estimate spectra of our z = 6.6 and 5.7 LAEs (thick red 
and black lines, respectively, in Figure fT7|) . the slope of 
red wing is found to be slightly sharper for z = 5.7 LAEs 
than z = 6.6 LAEs. This would be a hint of flattening 
of Lya line profiles from z = 5.7 to 6.6. A future study 
with better spectroscopic data would address the issue 
of this slight evolutionary effect. 

Interestingly, our composite spectra indicate that Lya 
emission may not be simple smooth asymmetric lines, 
but with a substructure. Comparing our composite spec- 
tra with the curves of instrumental spectral resolution in 
the bottom panels of Figure \TE[ we find knees in the blue 
tail of the composite LAE spectra at ~ 9185 and 8152A 
for z = 6.6 and 5.7, respectively, both of which corre- 
spond to the same rest-frame wavelength about ~ 0.7A 
bluer than their line peaks. Since similar knees are not 
seen in our sky line spectra (e.g. the bottom right panel 
of Figure [3]), these knees are not made by DEIMOS's 
instrumental line profile. Another possibility may be a 
residual of sky subtraction, but there are no clear reasons 
why the residual is mostly positively scattered in statis- 
tical sense to make the knees. The peaks of the knees 
are detected at the 3.5 and 2.2 sigma levels for LAEs 
at z = 6.6 and 5.7, respectively, at the same rest-frame 
wavelength, which is difficult to be explained by the ran- 
dom errors (see la errors associated with the composite 
spectra in Figure [T5|) . Thus, the knees in the blue wings 
of z = 6.6 and 5.7 spectra are probably real. 

5.2. Lya Velocity-Luminosity Relation 

Figure [T|)] presents a weak, but a positive correlation 
between AVfwhm and logL(Lya) in our z = 6.6 LAEs. 
A linear function fit to the z = 6.6 AVFWHM-logL(Lya) 
data indicates that the positive correlation is found at 
the 2.5(7 level. This tren d is opposite to the one fo und by 
iKashikawa et al.l ([20061 ) . IKashikawa et all ([20061 ) study 
LAEs with logi(Lya) ~ 42.3 — 43.0 that are fainter than 
ours by a factor of 2-4 on average. A distribution of data 
points of IKashikawa et al.l (|2006l ) is similar to ours in 
the luminosity range where both studie s have measure- 
ments. Figure 11 of IKashikawa et all (J2006I ) presents 
two faint LAEs that have large FWHMs, which ap- 
parently make the anti-correlation. If these two faint 
LAEs are largely up-scattered by statistical errors or 
sample va riance, the anti-co r relatio n is not clearly found. 
Moreover, IKashikawa et al.l (|2006l ) consider no velocity- 
width measurement errors in their linear-function fitting, 
while the uncertainties of velocity-width measurements 
increase towards faint luminosity. Thus, the previous 
conclusion of the anti-correlation is probably not strong. 
On the other hand, in Figure I16[ our positive correla- 
tion is apparent when our data of three bright LAEs in 
logL(Lya) ~ 43.3 — 43.6 are included. When these three 
data points are not used for our fitting, neither of posi- 
tive nor negative correlation is identified beyond the la 
level. If these three bright LAEs are not typical ones due 
to a sample variance, the clear positive correlation would 
not be found. Although there is a possibility that these 
three LAEs are not typical due to sample variance, these 
bright LAEs have small errors that are difficult to pro- 
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TABLE 6 
Average Lyo Line Properties 




<*> 

(1) 


(-^Lya) 

(10 43 erg s- 1 ) 
(2) 


AVfwhm (direct) 

(km s^ 1 ) 

(3) 


A VpwHM (Gaussian) 

(km s" 1 ) 

(4) 


Error(AVp WH M) 

(km s" 1 ) 

(5) 


6.56 
5.71 


1.39 + 0.06 
0.96 ± 0.03 


251 

260 




270 
265 


16 

37 



Note. — (1): Average redshift and (2): average Lya luminosity of our spectroscopic 
data. (3): Average FWHM velocity width given by the direct measurement, (4): average 
FWHM velocity width determined by the Gaussian fitting, and (5): error of FWHM velocity 
width measurement. 



duce largely up-scattered measurements. We conclude 
that there is no anti-correlation between Lya luminosity 
and line width at z = 6.6, and that, if our spectroscopic 
sample is not a biased one, Lya velocity width positively 
correlates with Lya luminosity at z = 6.6 in the lumi- 
nosity range of logL(Lya) ~ 42.6 — 43.6. 

6. DISCUSSION 

6.1. Constraints on the Cosmic Reionization History 

In this subsection, we discuss implications for cosmic 
reionization based on our LF, clustering, and Lya line 
profile results with the aid of theoretical models. Be- 
cause some conclusions of theoretical models would de- 
pend on their assumptions and methods such as analyti- 
cal, semi-analytical, and numerical techniques including 
radiative transfer, we compare our observational mea- 
surements with as many various theoretical models avail- 
able to date as possible. In this way, we aim to obtain 
implications for cosmic reionization with less model de- 
pendencies. 

6.1.1. Evolution of Lya LF and Luminosity Density 



In Section I3TI we have found that Lya LF decreases 
from z = 5.7 to 6.6 at the > 90% confidence level, and 
that the decrease is 30% in luminosity for the case of 
pure luminosity evolution. Because Lya LF evolution is 
made not only by cosmic reionization but also by galaxy 
evolution, we should interpret the decrease of Lya LF 
carefully. In fact, UV-continuum luminosity function of 
dropout galaxies also decreases from z = 6 to 7 — 8 in ob- 
servat i onal data (e.g. lOuchi et al.ll2009bt IBouwens et al.1 
l2010at ICastellano et~aTll2010r i. Some fraction of Lya LF 
decrease may be explained by this galaxy evolution ef- 
fect. Figure [18] presents a number density ratio of z = 6.6 
to 5.7 LAEs as a function of Lya luminosity. We find 
no significant dependence on luminosity for the ratio 
within errors, i.e. by a factor of ~ 2 at the luminos- 
ity range of logi = 42.5 — 43.5. However, it appears 
that the ratio is relatively smaller at the bright lumi- 
nosity (logL = 43.0 — 43.5) than the faint luminosity 
(log-L = 42.5-43.0). Note that Figure H8l misses a mea- 
surement at the bin of logL = 43.3, because no LAEs 
are found at z = 6.6 in this bin. If the trend of relatively 
small number density of z = 6.6 LAEs at the bright 
luminosity is true, galaxy evolution may be more domi- 
nant on the LF evolution between z = 5.7 and 6.6 than 
the cosmic reionization effect. Since brighter LAEs are 
likely hosted by more massive dark halos whose forma- 
tion epoch is near these redshifts (see halo mass func- 
tions, e.g., ISheth fe Tormenl I1999D . evolution of num- 
ber density of bright LAEs would be more affected by 



galaxy evolution. On the other hand, simple models of 
LAEs suggest that brighter LAEs have a less reduction of 
Lya fl ux due to large ionized bubbles surrounding bright 
LAEs (|Haimanll2002l ). which would imply that evolution 
of bright LAEs are milder than that of faint LAEs in the 
reionization-effect dominant case. 

We need to quantify how much decrease of Lya LF is 
contributed from cosmic reionization or galaxy evolution. 
If we have an assumption of no galaxy formation effects 
in Ly a LF evolution, which is made in previous studie s 
(e.g. iMalhotra fc Rhoaddl200l iKashikawa et al.ll2006| ). 
a ratio of IGM's Lya transmission at z = 6.6 to the 
one at z = 5.7 is 0.7. However, it is not clear whether no 
galaxy evolution assumption is correct. UV LF evolution 
of LAEs between z = 5.7 and 6.6 could resolve the degen- 
eracy of Lya LF evolution bet ween cosmic reion i zation 
and galaxy formation. Although IKashikawa et all (J2006I ) 
claim that there is no evolution of UV LF of LAEs be- 
tween z = 5.7 and 6.6 based on their Lya-subtracted z' 
photometry, we have found that UV LF of LAEs cannot 
be derived from Lya-subtracted z' photometry with an 
accuracy better than a factor of 2-3, due to large uncer- 
tainties. In fact, UV magnitudes estimated with the Lya 
subtracted z' photometry have large errors, > 0.4 — 0.7, 
as shown in Table [2] These large errors arc raised, be- 
cause fluxes in z' band are dominated by the strong Lya 
line, but not by the faint UV continuum. Only with opti- 
cal z' band photometry, UV LFs of LAEs can be reliably 
derived up to z < 6 LAEs whose Lya lines do not enter 
the z' band (e.g-jHu et alj|2004 IShimasaku et al.|[200l 
lOuchi et al.l 120081 ) . and near- infrared photometry is re- 
quired to derive reliable UV LF of LAEs at z = 6.6. It 
should be noted that, by the same reason, EWq values 
of UV-continuum faint LAEs (EWq > 100 A) are very 
poorly constrained (Table [2]), and that a number- EWq 
distribution of our z = 6.6 LAEs is not obtained reliably. 
Since we cannot use an uncertain UV LF of LAEs at 
z = 6.6, we investigate effects of galaxy formation with 
the other methods. 

We calculate luminosity densities, p, of Lya lines 
from LAEs and UV continua from dropout galaxies, and 
present ratios of p z / Pz~e in the top panel of Figure [T9l 
where p z and p z ^e ar e luminosity densities of redshifts z 
and 6, respectively. For these calcul ations, we use Lya 
LFs at z = 6.6 (this study) and 3. 1 (lOuchi et all 1200 
as well as UV LFs at z = 2 - 7 (IBouwens et all 1200 .. 
iReddv fc SteideH 120091 lOesch et al.ll2010t ). We estimate 
two sets of p from integrals down to the observed lumi- 
nosity (log L Ly a = 42.4 for Lya or Mrjv = -18 for UV) 
and down to zero luminosity; the latter is probably near 
a total p. We confirm that the ratios of these different 



20 



Ouchi et al. 



i.O 



II < 



to a 
II < 



0.5 



II » 
-II 

1 II II 



42.5 43.0 43.5 

logZ(Lyoc) [ergs s" 1 ] 



Fig. 18. — Ratio of number densities of z = 6.6 LAEs to z = 5.7 
LAEs, as a function of Lya luminosity. Squares plot observational 
resul ts from LAEs at z = 6.6 (this study) and z = 5.7 i Ouchi et alj 
120081) . The data binning is the same as that in Figures 16171 

estimates agree within the error bars in Figure 1191 We, 
thus, refer the latter estimate for our fiducial results in- 
cluding no systematic bias from observations. The top 
panel of Figure H9l indicates that Lya p z /p z ~6 stays con- 
stant within the errors between z = 3.1 and 5.7, but 
there is a drop from z = 5.7 to 6.6 beyond the error bar. 
This drop of Lya p z / p z ~e is originated from the decrease 
of Lya LF in this redshift range. On the other hand, a 
ratio of UV p monotonically decreases from z ~ 3 to 7. 
The decrease of UV p z j ' p z ^§ from z ~ 6 to 7 suggests 
that the cosmic star-formation rate density (SFRD) of 
galaxies decline at this redshift range, and this cosmic 
SFRD decline would contribute to the decrease of Lya 
p z /p z ~6 from z = 5.7 to 6.6. 

We evaluate the effect of cosmic SFRD decline on the 
basis of evolution of UV luminosity density. We assume 
that Lya luminosity density, p Lya , is proportional to UV 
luminosity density, p uv , 



P 



Lya 



K T, 



Lya fhya UV 



IGM 



/; 



(21) 



where Tjqm is a transmission fraction of Lya through 
IGM and /iff" is a Lya escape fraction within a galaxy 
through their inter-ste llar medium (ISM) , /iff" depends 
on gas infall+outflow (ISantos et al1 F2004: Dii kstra et al.l 
I2007al iDiikstra fc Wvithell2010D. distribution of galactic 
hydrogen ([ .Zheng et al. 2009, 2010J), and dust obscuration 
(jDaval et al.ll2010h . k is a factor converting from UV to 
Lya luminosities, which depend on stellar population, 
i.e., IMF, age, and metallicity. Assuming that stellar 
population of LAEs is the same at z = 5.7 and a given 
redshift, z, we can obtain a ratio of T t q^ /iff" at z = 
5.7 to a redshift, 
densities; 
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redshift to the one at z ~ 6. Red circles and blue squares represent 
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integrated down to the observed luminosity (logLL ya = 42.4 for 
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ity. Open symbols are shifted by +0.1 along redshift for clarity. 
Red solid and dashed lines simply connect the Lya points of red 
filled and open circles, respectively. Note that the ratios at z ~ 6 
have no error bars, because this is a definition of the ratios. In- 
stead, the errors of luminosity density measurements at z ~ 6 are 
included in the ratios at the other redshifts (z ^ 6). Bottom panel 
plots ratios of total Lya escape fraction, F a / C a , as a function of 
redshift. The ratios are -F CS c a at each redshift divided by F ea ya at 
z = 5.7. Again, the open circles are those calculated with the lumi- 
nosity densities integrated down to the observed luminosity, while 
the filled circles are estimated from the total luminosity densities 
integrated down to zero luminosity. Magenta dotted lines present 
linear fits of F eS c° ! /-F ca ^"_ 5 7 and redshift to filled and open circles 



" esc.z 

in the low-rcdshift range of z 



3.1-5.7. 



where the indices, z and z = 5.7, show redshifts. Be- 
cause the product, T I( Jm /iff"; is a total of Lya escape 
fractions, i.e. a fraction of Lya photons escaping from 
IGM and ISM, we refer to this product as a total Lya 
escape fraction. We write 



p~Lya 



rpLya rLya 
Iimi Jcsc 



TGM J esc (23) 

for simplicity. With this definition, the left hand side 
of eq. d22j) is rewritten as F^/ C a z /F^ c a z=57 showing 
redshifts with the indices. The eq. ([22]) means that 

^osc|z/^osc < z=5.7 is determined by the ratio of Lya lu- 
minosity densities divided by the ratio of UV luminosity 
densities. Since the ratio of UV luminosity densities, 
pY^ /pi=5 7, plays a role of a correction factor of cos- 
mic SFRD evolution, we refer this to a cosmic SFRD 
correction factor. The bottom panel of Figure [19] plots 
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Fig. 20. — Lyo transmission through IGM at z = 6.6 that is normalized by the one at z = 5.7. Left panel shows our observational 



estimates of T, 



Lya 
IGM,z = 6.6 



nLyo, 



Lya\ 



/Tjq m z _ 5 7 . Black squares plot T : q m z _ g g/TjQ M z _ 5 7 in the case of no evolution of Lya escape fraction (/ CS c ) 



Black solid and dashed lines are the best estimate and ±1<t errors of T, 



Lya 



pLya 



rLya 



'lGM,z=6.6/ T lGM,z=5.7' respectively, for the no f eB y c evolution case. 
Similarly, gray triangles and line represent the lower limit in the case of the f e J£ a evolution. Right panel compares these estimates with 
model predictions. Gray region shows our observational constraints (ztlc) on T^^ z _„ g/Tj,?^ z _ 5 7 . Solid and dotted lines represent the 

model predictions of Diikstra ct al. (2007b) as a function of ionized bubble radius in comoving Mpc for the cases with and without ionizing 
background boosts of undetected surrounding sources, respectively. Ticks at top present the relation between t ypical bubble rad i us an d 
neutral hydrogen fraction of IGM for xjji — 0.4, 0.3, 0.2, and 0.1, which arc predicted by the analytic model of IFur lanctto ct al. (2006). 



Two sets of ticks in the right hand side of the panel denote the predicted relations between xjji and T t qj^ z=g 6 /Tj 
of lSantos et~aTI pOOl ) and lMcQuinn et al.l (2QH3). 



Lya 
IGM,z=5.7 



in the models 



FeZ%/ F lL%=5.7 as a function of redshift. We find that 
the total Lya escape fraction might show a possible de- 
crease by ~ 20% from z = 5.7 to 6.6, albeit with a large 
error. Thus, the evolution of F^ s ya z is small and less sig- 
nificant than that of the ratio of Lya p z shown in the 
top panel of Figure [19] On the other hand, F^ s ya z in- 
creases from z = 3.1 to 5.7 in the bottom panel of Figure 
[T9l This is probably because a fraction of Lya escap- 
ing from ISM (i.e. /^ Q ) increases from low-z towards 
z = 5.7. If one naively extrapolates the linear fit of the 
relation at z = 3.1 — 5.7 between Fl^/ F c y c a z=z 7 an d 
redshift towards high-z in the bottom panel of Figure 
[T9l F^ y c a at z = 6.6 is smaller than that of the extrap- 
olation (magenta dotted lines in Figure [T9| just beyond 



-,Lyas 



and 



the la error. However, the evolution of IGM {T lGM ) 
galaxies (fc/ c a ) cannot be clearly distinguished. 

Thus, we evaluate fl£. a , and estimate evolution 
of Lya transmission of IGM from z = 5.7 to 6.6 

with a quantity of T t ^, z =6.6/ T igm,z=5.7- Because 
the estimates based on p^ ya /p z ^ 7 shown above 
would provide additional uncertainties raised by the 
integration ranges for Lya luminosity densities, we 
perform calculations with ratios of Lya luminosity by 
replacing p\ ya j pJJ^ 7 with the ratio of Lya luminosity, 

L^ ya /L 2 ^g 7 , in eq. (|22|) at given data points of observed 
number density, although this improvement of calcu- 
lations provide only negligible changes in our results. 



Note that the quantity of T, 



Lya 



can be written 
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V esc, z— 6.6/ esc, z— 5.7/ w esc, z— 5.7/ J esc, z— 6. 



by 



the 



definition (eq. I23[) . Here, the unknown physical quantity 
is the evolution of the ratio of Lya escaping fractions 



of ISM, i.e., (/ CS c Q z=5 .7//esc Q z=6.6)- Wc consider the 
two cases of f^/ c a evolution under an assumption of 
no luminosity dependence of f^ a for simplicity. The 
first case is no evolution of f^ y c a between z = 5.7 and 



pLyct 

esc, z— 6.6/ 
Lya 
IGM,z=5.7 — 



= 1. We find that 
0.80 ± 0.18, and 



6.6, which gives (/ eS c a z=5 .7//< 
the ratio is Tj^ z=6 6 /T, 
that IGM transmission at z = 6.6 decreases at the la 
level. For the second case, we estimate it with a / C L S J° 
evolution between z = 5.7 and 6.6. We assume that 
the redshift evolution of /^ a between z = 5.7 and 6.6 
is the same as that of linear increase per redshift found 
in the low redshift regime of Z = 3.1 — 5.7. To quantify 
the evolution of /<J^ Q at the low-z (z = 3.1 — 5.7), 
we can rewrite the eq. 



rLya i r 

•/esc,z— 3.1/ J « 



Lya 
esc, z— 5.7 



(^ L sr Z =3.l/i ? e L srz=5.7)(r i L G y M, z=5 . 7 /T I L G y M,z= 3 .l)- Figure M 

shows that (^=3.1/^=5-7) = °- 48 - We estimate 
the ratio of (^^,2=5.7/^^,2=3.1) with tnG a verage 
GP optical depths. Assuming that IGM absorbs a 
blue half of symmetric Lya emission line, we obtain 
(^i L G y M,z=5.7/ ^ L GM,z=3.i) | = (0-53/0.85) based on the 
estimates of iFan et al.l (|2006l ). 15 Note that this 
number is a lower limit, because Lya lines are red- 
shifted from the systemic velocity by a few hundred 
km s" 1 and less absorbed by IGM dPettini et al.l 120011: 
iMcLinden et al.ll2010t ISteidel et al.ll2010l ). The real num- 

15 The estimate of T r( j^ differs only by 3 — 5% from those with 
IMadaul lfi995T ) and lMeiksint pOOl ) (see lOuchi et~alll200l ) . 
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Fig. 21. — Comparisons of our Lya LF at z = 6.6 with those given by theoretical models. Red circ l es and line are our best-estimate 
LF as shown in Figurc[7] Left panel shows LFs predicted with semi-analytic models of lKobavashi et al.l ||201CT ). Solid, dashed, and dotted 



lines are predicted LFs with T^J!^ z _ g q/T-^q^ z _ 5 7 = 1.0, 0.8, and 0.6, respectively. Right panel plots a LF reproduced by the radiative 



_IGM,z=6.6' IGM. 
transfer model of Iliev et al. (2008) with the applied shift of AlogL(Lya) 

bcr of (I^S „ 5 . r /IlX^8.i) is > (0-53/0.85), closer to 



+31.6. 



unity. Hence the value of / L S c Q z=3 .i//o L S c Q z=5.7 is > °- 30 
(= 0.48 x 0.53/0.85). Extrapolating the redshift evolu- 
tion of /^ Q to 6.6, we obtain the ratio of Lya escape 
fraction at z = 5.7 to 6.6, (/ e L S c^= 5 .7//ic^=6.6) > °- 8L 



Thus, T, 



Lya 



IT, 



Lya 



IGM,z=6.6/ - t IGM,z=5.7 



> 0.65 in this second case 



with the redshift evolution of f^. a . 

The left panel of Figure 
T' 



Lya 



presents 



6 e/^iGM z=5 7 f° r the nrs ^ an d second cases 
with black and gray lines, respectively. Figure [20] also 
plots the ratio as a function of Lya luminosity for each 
observational data point with squares and triangles. 
Note that the second case only gives the lower limit, 
and that the gray triangles in the left panel of Figure 
20l represent the lower limits. Figure [20] shows that 
these lower limits are consistent with the results of 
the first case within the — la error, and, thus, the 
second-case result is included in the first-case result of 



z ~ 3 — 6 is probably reliable. We plot their Lya LFs at 
z = 6.56 in three cases of T^ z=66 /T^ z=57 = 1.0, 
0.8, and 0.6 16 . Although the models under-prcdict 
number densities of LAEs at the faint end (logL < 42.8) 
for all the cases, shapes of the predicted LFs agree at 
the bright luminosity (log L > 42.8). Because their 
faint-end LF may be strongly affected by complicated 
star-formation and supernova f eedback processes de- 
pe nding on model assumption s (|Kobavashi et al.1 120101 : 
see iNagashima fc Yoshiiir2004l for their models of galaxy 
formation components), we compare only their bright- 
end LFs. The left panel of Figure [21] shows that their 
bright-end LF of the T^ z=66 /T^ z=57 = 0.8 model 
reproduces our observational results. This constraint on 



r, 



Lya 



/Tu 



IGM,z=6.6/ J IGM,z=5.7 



is consistent with our estimate 



riLya 
L IGM,z=6.6 



m 



Lya 
IGM,z=5.7 



0.80 ± 0.18. It indicates that 



the IGM transmission at 



6.6, T, 



Lya 



IGM,z=6.6' 



is smaller 



than the one at z = 5.7 by 20% but just beyond the 
la error. We conclude that there would exist a small 
decrease of IGM transmission at the la level probably 
contributed by cosmic reionization, and that the ratio of 



L 



Lya 



IGM,z 

We 



=6.6/^] 



Lya 
IGM,z=5.7 



is ~ 0.8 and no smaller than 0.6. 



obtain another independent constraint on 

T roM,z=6.6/ T iGM,z=5.7 h y comparison with physical 
models including galaxy evolution and cosmic reioniza- 
tion effects. Left panel of Figure [21] co mpares Lya LFs 
xedicted with semi-analytic m odels of iKobavashi et all 
20101) . Because the model of Ikobavashi et al.l (|2010l) 
reproduces observed Lya LFs between z = 3.1 and 
5.7, the galaxy evolution component of this model at 



of T iGM.z=6.6/ r iGM,z=5.7 = 0.80 ± 0.18. Right panel of 
Figure [21] presents Lya LF predicted with the radiative 
transfer model of llliev et all (j2008| ). Because their 
luminosity is arbitrary, we have applied a shift of 
AlogL = +31.6 to their model that is roughly matched 
to our observational data points. The LF shape of 
the model agrees with that of our observations in the 
luminosity range of logL ~ 42.5 — 43.5, although there 
is a hint of a steeper LF slope in the model than our 
observational measurements. 

Next, we place a constraint on a neutral hy- 
drogen fraction of IGM, xhi, and a typical ion- 
ized bubble radius at z = 6.6 based on our esti- 
mates of r^J z=6 6 /T IGM ,z=5.7 - 0.8. Right panel 
of Figure [20] ticks xhi values at the corresponding 



16 Because IKobavashi eT"aH <20Ul) define T^™ z=5 7 = 1.0, 
these three models correspond to their models of T^ ya = 1.0, 0.8, 
and 0.6. 
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iLya 



^igm z=6 6/-^iGM,z=5.7 obtained by analytic models of 



Santos et al.l (I2004T) and radiative transfer models of 
McQuinn et al.l (|2007|) . In these comparisons, we assume 



that Tigm is contributed by scattering of Lya damp- 
ing wing of hydro gen IGM alone. In the models of 
ISantos et al.1 (|2004h . we have applied the realistic mod- 
els with a Lya line redshifted by 360 km s _1 from a 
systemic redshift, since redshifted Lya lines by a few 
hundre d km s" 1 are observationally found n ot only in 
LBGs (jPettiniet all [200 1; Steidel et all 12010ft. but also 
in LA Es (jMcLinden et al.ll2010h . In the IMcQuinn et al.l 
(|2007t) models, we use ratios of intrinsic to observed Lya 
luminosities at the cumulative number density of ~ 10~ 5 
Mpc -3 , which corresponds to logL ~ 43 near to the lu- 
minosity range of our observations. Both of the models 

0.8 corresponds to 



^5.7 



indicate that T^^ z=66 /Ti GM , z -- 
cehi — 0.2. Even with the errors of T, 



IGM,z=6.6/ J IGM,z=5.7, 



sAl 



we find that xjji is smaller than < 0.4 — 0.5. 

The right panel of Figure [JO] also presents cos- 
mic reionization mode ls of IDiikstra et al.1 l)2007bl) and 
IFurlanetto et al.l (|2006t l. The models of |D]jkstra_eLa]J 
(|2007al ) imply that our observational estimate of 



T, 



Lya 



6 6/^iGM z=5 7 — 0-8 corresponds to a typical ion- 
ized bubble s ize of > 40 comoving Mpc. Following the 
discussion of IDiikstra et al.l (12007IJ). w e apply the an- 
alytic model of IFurlanetto et al.l (|2006[ ) to this typical 
ionized bubble size. At the top of the right panel of 
Figure 1201 we tick the volume averag ed neutral fraction 
predicate d bv IFurlanetto et al.1 (120061 ). We find that the 
model of IFurlanetto et al.l ( 20061 ) infers that our lower 
limit of > 40 comoving Mpc corresponds to the neutral 
fraction of x^ < 0.1. 

In summary, all of these various theoretical mod- 
els, i.e. analytic, semi-analytic, and numerical 
models, indicate that our observational estimate of 



r, 



Lya 



iGM,z=6.6/ r iGM,z=5.7 = °- 80 ± °- 18 corresponds to 
%Ht 5, 0-2 ± 0.2. In the p revi ous studies of Ly a LFs, 
iMalhotra fc Rhoadsl (|2004l ) and iKashikawa et all (|2006l ) 
place upper limits of neutral fraction for xhi < 0.3 and 
< 0.45, respectively, with the comb ination of the i r mea - 
surements and one analytic model of lSantos et al.l (|2004f ). 
Our result is consistent with these previous studies, but 
our constraint is stronger and more robust than these 
previous results because of our better LF evolution de- 
termination with cosmic variance errors and the inclu- 
sion of model dependent errors based on various indepen- 
dent reionization models. If all of our deep fields are not 
strongly biased by patchy IGM distribution by chance, 
our conclusion of IGM neutral fraction, a;^ ;$ 0.2 ± 0.2, 
is not significantly changed. 

6.1.2. Evolution of Clustering 

Theoretical models predict that a clustering ampli- 
tude of observed LAEs is boosted, due to the addi- 
tional clustering of LAEs whose Lya photons can es- 
cape in the case that these LA Es reside in an ionized 
bubble at the reionization epoch (IFurlanetto et al.ll2006l: 
McQuinn et~aTI 120071: iLidz et al.lT2009l; cf. Illiev et al.l 
2008) . In Figure Q21 we find no sudden rise of bias from 
z = 5.7 to 6.6 beyond the error bars. This indicates that 
clustering of z = 6.6 LAEs is weakly affected by cos- 
mic reionization. With this result on bias evolution, we 
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Fig. 22. — Comparison with our clustering measurem ents and 
reionization models. This plot is the same as Figure 1121 but 
with the m odel predictions of an gular correlation function in the 
top panel { McQuinn et al. 2007) and bias in the bottom panel 
IIFurlane tto^et^dJ [2006J) . Top : Solid lines are the model predic- 
tions of McQuinn et al. (20Q3) for z = 6.6 LAEs whose hosting 
dark halo mass is 3 X 10 1(1 Mq. Neutral fractions of the universe 
are xm = 0.0, 0.3, 0.5, and 0.8 from bottom to top lines. Dotted 
lines indicate z = 6.6 LAEs with a dark halo mass of 7 X IQ 10 M@ 
(bottom) and 1 X 10 Mq (top) in the c ase of xhi = 0.0. B o ttom : 
At the right side of this plot, we present [Furlanetto e t all s 1120061 ) 
predicted bias values under the assumptions described in the text. 
For clarity, we also plot the best-bias estimate and la ranges for 
our bright (7VZ?921 < 25.5) LAE subsample with gray solid and 
dashed lines, respectively. 



constrain cosmic reionization models of LAE clustering 
predictions. 
Figure |2"21 compares the ACFs and bias of our LAEs at 
6.6 wit h theoretical pred i ctions of 



d. 
d. 



McQuinn et al 
2007| and IFurlanetto et all (|2006t) . IMcQuinn e~ 
2007f ) models are based on radiative transfer simula- 
tions. In their models, Lya emission fluxes are assigned 
to dark halos with a mass above a minimum halo mass. 
We find that the average hosting dark halo masses of 
LAEs are about 10 11±:l Mq over z = 2 - 7 in Section |431 
This result is also correct at only z = 2 — 6, even if we 
omit the results of z = 6.6 LAEs whose clustering might 
be contaminated by the cosmic reionization effect. If we 
assume that hosting halo masses of LAEs do not change 
significantly from z ~ 6 to 6.6 so as true for z = 3 — 6 
LAEs, the possible minimum halo mass of z = 6.6 LAEs 
is about 10 10 Mp,. The top pan el of Figure [25] plots the 
models of lMcQuinn et all (|2007f ) with the minimum halo 
mass of 3 x 10 10 M Q which is the closest to 10 10 M Q and 
the smallest halo mass in their models. The compar- 
ison of our z = 6.6 LAEs and these models indicates 
that the neutral hydrogen fraction of IGM is less than 



xni < 0.5 at 



6.6. We also plot the rest of models 



with masses of 7 x lCr " Mp, and 1 x 10 n M Q provided by 
IMcQuinn et al.l (|2007|) . and confirm that these massive 
halo models generally give a strong clustering, and that 
their constraints on xhi are all consistent with this rela- 
tively weak upper limit of Ihi ^ 0.5. Thus, the compar- 
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ison with iMcQuinn et al.l (|2007l ) model s places an upper 
limit o f xhi ^ 0.5. The comparison with lFurlanetto et al.l 
( 2006() models is shown in the bottom panel of Figure 
[22] We use r atios of small-scal e bias to hosting halo bias 
for LA Es in IFurlanetto et al.l (|2006t ). IFurlanetto et al.l 
(|2006| ) define the small-scale bias as a bias at a scale 
smaller than a typical ionized bubble size at the red- 
shift. At z = 6.6, near the end of reionization, a typical 
size of ionized bubbles is > 10 — 100 Mpc. Because our 
bias estimates of z — 6.6 LAEs arc mostly made at the 
scale of ~ 1 — 10 Mpc (see Figure IT2J) . we, thus, regard 
our bias val ues as a small-scale bia s that is defined in 
the study of IFurlanetto et all ([20061 ). For a hosting halo 
bias, we assume again no evolution of halo mass, and use 
b = 3.7 that is a bias of 10 10 M Q at z = 6.6 given by 
cq. (|19[) . Comparing the la limits of average bias for 
all sample (black dashed line) and the bright subsample 
(dark dashed line) in the bottom panel of Figure [22] we 
find that a neutral hydrogen fraction of IGM should be 
smalle r th an xhi 5, 0-5 at z = 6.6. Either IMcQuinn et al.l 
(120071) or IFurlanetto et al.l (120061 ) models would imply 
that a neutral hydrogen fraction is xhi ^ 0.5 at z = 6.6 
based on our clustering estimates. 

6.1.3. Evolution of Lya Line Profile 

Lya damping wing of neutral hydrogen IGM absorbs 
a Lya emission line at > 1216A in the partially and fully 
neutral universe. Because Lya damping wing absorption 
at > 12 16 A monotonically weakens towards red wave- 
lengths, the damping wing absorption changes a shape 
of Lya line profile. Since a Lya flux absorbed by Lya 
forest has an asymmetric profile with its flux peak near 
1216A, i.e. in the bluest wavelength of high-z Lya line, 
the absorption component of Lya damping wing (hav- 
ing a stronger attenuation in bluer wavelengths) should 
bro aden the Lya l i ne. We use models of Lya lines made 
by iDiikstra et al.l (]2007al ). and evaluate the Lya line 
broadening. We find that an FWHM of Lya line broad- 
ens by 7 — 10% in a very neutral universe where a LAE 
resides in an ionized bubble with a radius of 0.7 physi- 
cal Mpc (5 comoving Mpc at z = 6.6), compared with 
a LAE in more ionized universe with an ionized bubble 
radius of 2 — 10 physical Mpc (15 — 80 comoving Mpc at 
Z = 6.6). Because the curvature of Lya damping wing 
absorption is not steep, the Lya line broadening is not 
large, only by a few-10 percent level. Thus, high quality 
spectra arc needed to identify the Lya line broadening, 
if any, at an observed redshift. In §5.1| we obtain the 
direct (Gaussian) FWHM velocity widths of 251 (270) 
±16 km s" 1 at z = 6.6 and 260 (265) ±37 km s" 1 at 
z = 5.7. These measurements include errors of 6% and 
14%, which allow us to marginally identify a few-10 per- 
cent level broadening of Lya. Our data show that there 
is no large evolution of Lya FWHM beyond ~ 14%, but 
it is not clear whether a few-percent level FWHM broad- 
ening exists. In fact, there is a possible hint of flattening 
of Lya profile from z = 5.7 to 6.6 f £J5.1[) , although its 
significance level is below the la level. One of the goals 
of future LAE studies will be testing the broadening of 
L ya FWHM dow n to a few percent level. 

lHaiman fc Cenl (|2005| ) suggest a test of reionization 
with a relation between Lya luminosity and line width. 
They claim that LAEs in a neutral universe have an 
anti-correlation between Lya luminosity and line width, 



because fainter LAEs preferably residing in smaller HII 
regions are affected by a stronger damping wing at- 
tenuation than brighter LAEs. We c ompare our ob- 
servat ional results with the model of lHaiman fc Cenl 
(|2005[ ) in the luminosity range and velo city width sim- 
ilar to our observations. The model of lHaiman fc Cenl 
(|2005l) indicates that, if the universe is fully neutral 
(xhi = 1), an FWHM velocity width (line width) de- 
creases only by 30 km s _1 (lA) from logL(Lya) = 42.5 
to logL(Lya) = 43.5 for LAEs with a velocity width 
of 300 km s . In their model, no anti-correlation, 
but rather positive-correlation, is found for more ion- 
ized universe with a neutral fraction of .thi = 0.5 and 
0.25 in this relatively bright luminosity range. Thus, 
in this luminosity range, they predict that a weak anti- 
correlation appears only when the universe is nearly neu- 
tral. Our observational results indicate that there is no 
anti-correlation between Lya luminosity and FWHM ve- 
locity width ( §5.2p . Our best- fit function presented in 
Figure [16] shows a rather positive correlation and an in- 
crease of FWHM velocity width by 111 ±45 km s _1 from 
log L (Lya) = 42.5 to log L (Lya) = 43.5. If the model 
prediction of lHaiman fc Cenl (|2005| ) is correct, our ob- 
servations rule out the fully neutral universe at z = 6.6 
at the ~ 3cr level. 

In §5.11 we find an interesting knee feature in Lya line 
profiles of LAEs at z = 6.6 as well asz = 5.7. It would be 
suggestive of galaxy outflow or proximity effect, but it is 
not clear whether it can be interpreted with a reasonable 
physical picture. Because we find these knee features not 
only at z = 6.6, but also z = 5.7 when the universe is 
highly ionized (xhi ~ 10~ 4 ), the knee features would not 
be related to cosmic reionization but galaxy formation. 
The knee features may be important for understanding 
dynamics, UV radiation field, and structure of LAEs at 
high redshifts. 

6.1.4. Cosmic Reionization 



In Sections I6.1.1ll6~l~3] we have obtained the con- 
straints of neutral hydrogen fraction of xhi < 0.2 ± 0.2 
and iehi ^ 0.5 at z = 6.6 from the evolution of Lya 
LF and clustering, respectively, and ruled out the fully 
neutral universe at Z = 6.6 at the ~ 3er level by Lya 
line profiles. It should be noted that these three con- 
straints agree, even though these results are given by 
three independent observational quantities, i.e. Lya LF, 
clustering, and Lya line properties. In Figure I23[ wc 
plot our two relatively strong constraints of Xhi from 
Lya LF and clustering, and compare with the previous 
estimates and theoretical models. Our results are con- 
sistent with the previous ones based on LAEs as well as 
GRBs 17 and QSOs, and our constraint from Lya LF, 
^hi ;$ 0.2 ± 0.2, is generally stronger than the previous 
ones at z ~ 6.5 — 6.6. This upper limit at z = 6.6 prefers 
an early (z > 7) reionization on average, which requires 
more ionizing photons at z > 7. 

Triangle in Figure [23] is the la lower limit of red- 

17 The estimates of a single GRB would include an additional 
syste matic error of <5xhi ~ 0.3 owing to the patchiness of reioniza- 
tion (McQuinn et al. 2008). Since the z ~ 6.7 GRB estimate of 
IGreiner et al.1 ||2009H allows x m = 0.001 - 1 within the 1 .2<r level 
due to the strong degeneracy betwee n absorptions of a d amped 
Lya absorber and IG M, the result of Grcincr et al. (2009) is not 
included in Figure [23] 
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shift given by the WMAP7 result i n the instantaneou s 
reionazation case (z = 10.5 ± 1.2; lLarson et a.1.1 l201Qf ). 
This rules out the instantaneous reionization at z < 9.3, 
i.e. 100%-neutral fraction at z < 9.3 , at the la level. 
Figure [23] sho ws predictions of semi-a nalytic reioniza- 
tion models of IChoudhurv et al.l ()2008| ). For sources of 
reionization, they assume three minimum halo masses 
for star-forming galaxies, ~ 5 x 10 5 M Q , ~ 10 8 M Q , and 
~ lO 9 A/ at z = 6 for mini, smal l, and large halo 
cases, respectively. IChoudhurv et all (|2008| ) have con- 
cluded that reionization only by large halos is ruled 
out by the constraints of the electron sca ttering op- 
tical depth (WMAP3; ISpergel et all I2007J ) and vari- 
ous neutral hydrogen fraction upper limits available at 
the time of their study. Our results support their 
claim, and our relatively strong upper limit of Xjn ^5 
0.2 ± 0.2 favor their model including minihalos. Mini- 
halos can accomplish the relatively early reionization, 
as suggested by various theoretical studies. On the 
other hand, a change of star-forming galaxy properties 
is also possible. These could be an increase of escape 
fraction of ionizing p hotons such sugges t ed by high- z 
dropout observations (fOuchi et al.1 l2009bl: iBunker et al.1 
I2009t iFinkelstein et al.l 120091: iBouwens et alj|2010b[ ). It 
is also possible that t he changes of IMF and metallic- 
ity are i mportant (e .g. IStiavelli et al.l 120041 ). In fact, the 
model of lCenl (J2003I ) shown Figu re 123] i nclud es the effects 
of IMF and metallicity changes. ICenl (|2003|) argues that 
reionization took place twice; the first at z ~ 15 — 16 
by metal-free Population III stars with a top-heavy IMF 
and the second at z ~ 6 by Population II stars, which 
give a medium large Thomson scattering optical depth 
of 0.1 ± 0.03 similar to WMAP7 results (jLarson et al.l 
12010( 1. This scenario gives a fairly low neutral fraction 
at z = 6.6 that is consistent with our LAE constraints. 
Although it is not clear whether reionization took place 
twice, metal-free star-formation would characterize the 
history of cosmic reionization. 

6.2. Role of LAEs in Galaxy Formation History 

The evolution of LAE LF, clustering, and Lya line 
profiles includes various hints for understanding LAEs 
in galaxy formation history. As shown in Figure 1141 a 
typical bias of LAEs is smaller than that of LBGs over 
z = 2 — 7. Thus, a typical halo mass of LBG is estimated 
to be ~ 1O 12±1 M0, a bout one order of magnitude larger 
than that of LAEs (iHamana et all 120041 lOuchi et all 



2004b, 20 0~5bl; iLee et al.ll2006L 120091: iMcLure et al.ll2009t 
Hildebrandt et aLll20(M) . 18 It is possible that some 



galaxies become LAEs with a halo mass of ~ 10 ll M & 
then evolve into LBGs with a halo mass of ~ 10 Mq 
via mass assembly such as mergers and accretion. If this 
scenario is true, a typical LAE could be a progenitor of 
typical LBG. In this scenario, two observational ques- 
tions are consistently answered. The first question is the 
different trend of LBG and LAE evoluti on. The LF of 
LBGs decreases from z = 3 to 6 (e.g. IBouwens et al.l 

18 Note that LBGs referred here are typical dropout galax- 
ies with a UV luminosity near L* , and that, in our discussion, 
we do not include the population of the LAE analog of dropout 
galaxies with a ve ry faint UV magnitude such those studied by 
IStark et al.l 1 120101) . Similarly, LAEs are typical LAEs so far ob- 
served in logL(Lya) ~ 42 — 44, and not those fainter than this 
luminosity range. 
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Fig. 23. — Neutral hydrogen fraction, ajjji, of IGM as a 
function of redshift. Top and bottom panels are the same, but 
with a vertical axis of linear and log scales, respectively. Filled 
square and circle are the upper limits of xjji that we obtain from 
the evolution of Lycr LF and clustering, respectively. Open di- 
amond and pen tagon denote the upper li mits from Lyq LF at 
z = 6 .5 given by [M alhotra & Rhoads (2004) and Kashikawa ct al. 
1120061 ). Open hexagon is the upper limit estimated f rom the con- 
strain ts of Lyo damping wing of GRB at 2 = 6.3 IjTotani et al.l 
120061 ). Filled hexagon and pentagons indicate constraints given by 
GRB spectra (Gallcrani ct al. 2008b) and QSO dark gap statistics 
IjGallera ni ct al. 2008a), respectively. Filled diamonds represent 
the m easurements from GP optical depth of SDSS QSOs dFan et al.l 
120061 ). Triangle plots the la lower-limit of re dshift of a neutral 
universe given by WMAP7 Larson ct al. (2010) in the case of in- 
stantaneous reionization. Avoiding overlapping symbols, we give a 
small offset along redshift to the positions of the filled circle and 
the open diamond. Dotted, dashed, and solid lines show the evolu- 
tion of xjji f or minihalo, small , and large halo cases, respectively, 
predicted by Choudhury ct al. (2008). Gray solid line presents the 

redic tion in the double reionization scenario suggested by ICenl 

20031 ). 



|2007| ). while the LF of LAEs does n ot change much in 
z = 3 — 6 (e.g. lOuchi et al.l [20 08). In this scenario, 
LAEs would form earlier than LBGs, and the signifi- 
cant drop of LF may be found only for LBGs in this 
intermediate redshift of z = 3 — 6. This is because the 
number density of LBG's massive (~ 10 12 M Q ) hosting 
halos sharply drops from z = 3 to 6, while that of LAEs' 
less-massive (~ 1Q 11 Mq) hostin g halos does not decreas e 
much based on the halo model (Sheth & Tormcn 1999). 
The redshift of z = 3 — 6 might be still the major forma- 
tion epoch of LBGs, while the major formation epoch of 
LAEs would be earlier than LBGs, i.e. at z > 6, near 
the epoch when we find a decrease of LAE LF (from 
z = 5.7 to 6.6) in this study. The second question is 
the deficit of strong Lya emitting galaxies among a UV 
bright population. It is found that UV -bright galaxies do 
not have a strong Lya emission line (|Ouchi et alj 12003 : 
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Ando et all 120061: IShimasaku et all 120061: IVanzella et al.l 
20071120091 fOuchi et al.ll2008t l. This trend is usually dis- 



cussed on a plane of UV luminosity and Lya EW, and 
no objects ar e found in the UV -luminous and EW-large 
regime (e.g. I Ando et"aT1l2006l ). In the case of our sce- 
nario, a typical LAE evolves into a typical LBG, i.e. less- 
massive galaxies with a strong Lya emission line become 
massive galaxies with a weak or no Lya emission line. 
When this evolution is investigated on this plane, objects 
move from UV-faint to UV-luminous regimes, because 
UV luminosity positively correlates with stellar mass at 
high redshi fts jPapovich et al.ll200ll : iSawicki et a7]|2007t 
lYabe et all 12009ft . The UV-faint (less-massive) objects 
are, first, distributed widely in EW values on this plane 
before the evolution. Then, these objects end up in 
the UV-bright and small-EW regime on this plane af- 
ter the evolution. The tendency of Lya deficit among 
UV-bright galaxies would be consistent with the sce- 
nario of this LAE-LBG evolution sequence. In this way, 
this evolutionary scenario would provide the answers to 
these two observatio nal questions. The recent study of 
IVanzella et al.l (|2009h has found that LAEs have a more 
compact UV-continuum morphology than LBGs with no 
Lya emission at z ~ 4. If our scenario of LAE-LBG 
evolution is correct, star-formation activities start at the 
center then extend to the outskirts of galaxies, which is 
suggestive of the inside-out picture of galaxy formation. 
The dotted lines in Figure ITBI show evolutionary tracks 
of dark halos for galaxy-conserving models. The galaxy- 
conserving model assumes that the motion of galaxies is 
purely driven by gravity, and that merging does not take 
place. In this case, the bias value of galaxies decreases 
as the Universe evolves with time, 

b s = l + (b° g -l)/D(z), (24) 

where &° is a bias at z = (|Fry||1996() . Under the as- 
sumption of the galaxy-conserving evolution, dark ha- 
los of z = 3 LAEs evolve into dark halos with b g ~ 
1.1 — 1.3 at z = 0. Because, i n more realistic exte nded 
Press-Shechter formalism (e.g. lLacev fc Cold 119931 ) . av- 
erage evolutionary tracks come slightly below those of 
the galaxy-conservi ng models (see, e.g., Figure 13 of 
llchikawa et al.ll2007|) , dark halos of z = 3 LAEs probably 
evolve into dark halos with b g ~ 1 which host present- 
day ~ L* galaxies in cluding Milky Way, as suggested by 
iGawiser etaTI (|2007f ). On the other hand, in Figure [131 
LAEs at z ~ 4 — 7 have a significantly larger bias than 
z — 3 LAEs, and the increase of the LAE bias measure- 
ments toward high-z (nearly along solid lines) is steeper 
than those of galaxy-conserving evolution (dotted lines). 
The galaxy-conserving model indicates that dark halos 
of z ~ 4 — 7 LAEs evolve into those with b g ~ 1.5 — 2 
at z = whose value is higher than those of z = 3.1 
LAEs (b s ~ 1.1 — 1.3), which implies that descendants 
of z ~ 4 — 7 LAEs are different from those of z = 3 
LAEs. On average, the descendants of z ~ 4 — 7 LAEs 
would be more massive than those of z = 3 LAEs. Ma- 
jority of LAEs at z ~ 4 — 7 are probably not ancestors 
of Milky Way, but today's large galaxies more massive 
than Milky Way, although there should exist some LAEs 
that become today's L* galaxies by a probability process 
of halo mass build up. These implications from our clus- 
tering resul ts are similar to those from theoretical pre- 
dictions by ISalvadori et al.l (|2010l) who claim that only 



~ 2% of Milky Way progenitors can be LAEs. 

7. CONCLUSIONS 

We have identified 207 LAEs at z = 6.6 in the 1 
deg 2 area of SXDS field down to L > 2.5 x 10 42 erg 
s~ 4 and EW > 14 A by deep and wide-field narrow- 
band and broad-band imaging of Subaru/Suprime-Cam. 
Nineteen Lya lines are confirmed by the high-quality 
Keck/DEIMOS spectra, and none of interlopers have 
been found in our sample through our extensive deep 
spectroscopy campaign as well as the SXDS project spec- 
troscopy consisting of 3,233 objects at z = — 6. We 
have obtained the Lya LF, ACF, and Lya line profiles 
to constrain cosmic reionization and early galaxy forma- 
tion with the aid of recent theoretical model predictions. 
We have also derived ACFs of LAEs at z = 3.1-5.7 
in the SXDS field to find the evolutionary trend of LAE 
clustering from z ~ 3 — 7 in the framework of the ACDM 
model. The major results of our study are summarized 
below. 

1. Our Lya LF of z = 6.6 LAEs shows the best-fit 
Schechter parameters of <j>* = 8.5+2i2 x 10 _4 Mpc~ 3 and 
Ll ya = 4.4±g;| x 10 42 erg s" 1 with a fixed a = -1.5, 
where the errors include uncertainties of statistics and 
cosmic variance. The combination of statistics and cos- 
mic variance errors presents scatters of number density 
measurements up to a factor of ~ 10 among 5 subficlds 
of ~ 0.2 deg 2 areas, although the typical scatters of the 
subfield number densities are not far beyond the errors 
of Poisson statistics. Comparing this LAE LF at z = 6.6 
with the one at z = 5.7, we find that the Lya LF de- 
creases from z = 5.7 to 6.6 at the > 90% confidence 
level. A more dominant decrease of L* (luminosity evo- 
lution) than cf>* (number evolution) is preferable. The 
decrease of Lya LF from z = 5.7 to 6.6 is ~ 30% in the 
case of pure luminosity evolution. Note that this 30%- 
luminosity decrease is too small to be identified by the 
previous studies, due to their large uncertainties from 
small statistics and cosmic variance. 

2. We have identified a significant angular-correlation 
signal for our z = 6.6 LAE sample. This is the detection 
of clustering signal for the most distant galaxies, to date. 
The correlation length and bias are r = 2 — 5 Hf 00 Mpc 
and bias of b = 3 — 6, respectively. There is no sudden 
boost of clustering amplitude given by cosmic reioniza- 
tion at z = 6.6. In the framework of ACDM models, the 
average hosting dark halo mass inferred from clustering 
is 10 10 - lO n M . The duty cycle of LAE population, a 
product of star-formation and Lya emitting duty cycles, 
is a few 0.1 to a few percent, roughly ~ 1%, although the 
constraint on the duty cycle involves large uncertainties 
because of the bias estimate errors and the fairly flat re- 
lation between number density and bias of halos from the 
ACDM model. 

3. Based on our high quality DEIMOS spectra, we 
have found that most of Lya emission lines present a 
clear asymmetric profile with the average FWHM veloc- 
ity width of 251 ± 16 km s _1 (direct measurement) at 
logL(Lya) ~ 42.6 - 43.6, and that the average FWHM 
velocity width does not largely evolve from z = 5.7 to 6.6 
beyond errors of our ~ 40 km s _1 , which is dominated by 
the error of our reference z = 5.7 LAE. There is no anti- 
correlation between Lya luminosity and velocity width. 
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If our spectroscopic sample is not biased, Lya velocity 
width positively correlates with Lya luminosity at z = 
6.6 in the luminosity range of logL(Lya) ~ 42.6 — 43.6. 
We identify a knee feature in a blue tail of Lya line in 
our composite spectra of LAEs at z = 6.6 as well as 5.7, 
which cannot be explained by statistical and systematic 
instrumental errors. These knee features, if true, would 
be important for understanding dynamics, UV radiation 
field, and structure of LAEs. 

4. We compare evolution of Lya LF, clustering, 
and Lya line profiles from our observations with vari- 
ous reionization models including analytic, semi-analytic, 
and radiative transfer models. Although there would ex- 
ist a small ~ 20%-decrease of IGM transmission from 
z = 5.7 to 6.6 due to cosmic reionization, the compar- 
isons of all models and observational quantities reach the 
same conclusion that hydrogen IGM is not highly neu- 
tral at z = 6.6. The upper limit of neutral fraction is 
xhi < 0.2 ±0.2 from our Lya LF evolution and xhi < 0.5 
from our clustering evolution between z = 5.7 and 6.6. 
A fully neutral universe, xhi = 1, at z = 6.6 is ruled 
out by no large evolution of Lya velocity width from 
z = 5.7 to 6.6 and by no anti-correlation between Lya 
luminosity and velocity width at z = 6.6. All of these 
reionization tests with Lya LFs, clustering, and line pro- 
files, agree that a neutral hydrogen fraction of IGM is 
not high at z = 6.6. Our strongest constraint from Lya 
LF, xhi ^$ 0.2 ± 0.2, implies that the major reionization 
process took place early, at z > 7. 

5. Calculating ACFs, bias, and hosting dark halos of 
low- 2 LAEs of SXDS at z = 3.1 — 5.7 in the same manner 
as those of our z = 6.6 LAEs, we find that hosting dark 
halo masses stay at the similar value of 10 11±1 M Q over 
z = 3 — 7. It implies that LAEs are galaxies at the evolu- 



tionary stage for all or some type of galaxies whose dark 
halos have reached a mass of ~ 10 11±1 M Q . Because a 
halo mass of typical (~ L*) LBGs is ~ 10 12±1 M Q , about 
one order of magnitude larger than that of LAEs, there 
is a possibility that some galaxies become LAEs with a 
halo mass of ~ 10 u Mq in the ~ 1% duty cycle, then 
evolve into typical LBGs with a halo mass of ~ 10 12 A/ Q 
via mass assembly such as mergers and accretion. This 
scenario consistently explains two observational results 
of LBGs and LAEs; the differen t LF evolutionary trend 
between LBGs and LAEs (e.g. lOuchi et af] l2008h . and 
the deficit of strong Lya emitt ers among UV bright pop- 
ulation (e.g. lAndo et al.ll2006f ). 
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